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ABSTRACT 


A  new  symmetrical  number  system  with  applications  in  parallel  signal  processing  is 
investigated.  The  "Robust  Symmetrical  Number  System"  (RSNS)  is  a  modular  system  in 
which  the  integer  values  within  each  modulus,  when  considered  together,  change  one  at  a 
time  at  the  next  position  (Gray  code  properties).  Although  the  observed  dynamic  range  of 
the  RSNS  is  somewhat  less  than  the  optimum  symmetrical  number  system,  the  Gray  code 
properties  make  it  particularly  attractive  for  folding  analog-to-digital  converters.  With  the 
RSNS,  the  encoding  errors  (due  to  comparator  thresholds  not  being  crossed 
simultaneously)  are  eliminated,  as  is  the  need  for  the  corresponding  interpolation  signal 
processing  (reduced  complexity).  Computer  generated  data  is  used  to  help  determine  the 
properties  of  the  RSNS.  These  properties  include  the  largest  dynamic  range  (number  of 
distinct  consecutive  vectors),  and  the  position  of  the  largest  dynamic  range  within  the 
system.  The  position  of  the  maximum  unambiguous  dynamic  range  is  also  quantified. 
Least  squares  analysis  of  2  and  3  moduli  systems  is  used  to  formulate  closed-form 
expressions  for  the  dynamic  range.  To  compare  the  advantages  of  the  RSNS  with 
previously  published  results,  the  transfer  function  of  a  3-channel  RSNS  folding  analog-to- 
digital  converter  architecture  (mi  =3,  nt2  =  4,  and  m3  =  5)  is  numerically  evaluated  using 
SPICE. 
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IX 


I.     INTRODUCTION 

A.         MODULAR  NUMBER  SYSTEMS 

This  study  addresses  a  new  modular  number  system.  Before  going  into  the  details 
of  the  new  system,  it  is  useful  to  explain  the  properties  of  the  some  of  the  other  modular 
number  systems  and  examine  their  applications  in  Signal  Processing.  One  of  the 
drawbacks  of  the  modular  signal  processing  is  the  introduction  of  large  errors  (encoding 
errors)  in  the  decoded  binary  or  decimal  values  when  only  small  errors  occur  within  each 
modulus.  The  encoding  errors  and  their  elimination  are  the  main  reason  for  defining  a 
new  modular  number  system. 

1.         Residue  Number  System 

The  residue  number  system  (RNS)  is  representative  of  a  sawtooth  folding 
waveform  with  folding  period  equal  to  the  modulus  m.  An  introduction  to  the  RNS  has 
been  given  by  Szabo  [1].  It  has  been  established  as  an  important  tool  in  parallel 
processing  applications  and  high-speed  computations.  The  RNS  can  serve  as  a  source  for 
parallel  operations  (moduli)  that  are  of  smaller  computational  complexity.  Each  parallel 
operation  has  a  different  modulus  and  requires  only  a  precision  in  accordance  with  that 
modulus.  The  folding  waveform  of  the  RNS  is  generated  as 

a  =  r  (mod  m)  ( 1 ) 

where  0  <  r  <  m,  where  m  is  the  modulus  (channel  value)  which  are  PRP  and  r  is  the 
residue  of  a  modulo  m.  Figure  1.1  shows  the  RNS  folding  waveforms  for  mi  -  4  and  m2 
=  5.  Going  from  integer  values  within  each  modulus  to  decimal  value  requires  the 
Chinese  Remainder  Theorem  (CRT)  [2-5].  The  integers  r  within  each  channel  (modulus) 
are  representative  of  a  folded  waveform  with  the  period 

Prns  =  rn  (2) 

The  dynamic  range  (region  of  no  ambiguities)  of  the  RNS  is 

N 
MRNS=Ylm,  (3) 
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and  m,  is  the  i    pairwise  relatively  prime  (PRP)  modulus  (channel  value),  1  <i  <N. 
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Figure  1.1:  The  Folding  Waveform  of  the  RNS  for  mi  =  4  and  m2  =  5 


2.         Symmetrical  Number  System 

The  SNS  is  also  composed  of  a  number  of  pairwise  relatively  prime  (PRP) 
moduli.  The  integers  within  each  SNS  modulus  however  are  derived  from  a 
symmetrically  folded  waveform  [6].  The  integer  values  within  each  SNS  modulus  are 
derived  from  a  symmetrical  folding  waveform  and  therefore  incongruent  modulus  m. 
Due  to  the  presence  of  ambiguities,  the  set  of  integers  within  each  SNS  modulus  does  not 
form  a  complete  residue  system  by  themselves.  It  is  well  known  that  the  inclusion  of 
additional  redundant  moduli  can  effectively  detect  and  correct  errors  that  may  occur 
within  a  RNS  representation  of  a  number.  The  SNS  formulation  is  based  on  a  similar 
concept,  which  allow  the  ambiguities  that  arise  within  the  SNS  to  be  resolved  by  using 
various  arrangements  of  the  SNS  moduli.  The  integer  values  within  each  SNS  folding 
waveform  is  given  as 
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0,1,..., 
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,  when  m  is  odd 


,  when  m  is  even 
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where  |_*J  indicates  the  greatest  integer  less  than  or  equal  to  x.  Both  have  size  1  x  m 
and  consist  of  the  symmetrical  remainder  elements  Xh,  0  <  h  <  m.  Figure  1 .2  shows  the 
SNS  folding  waveform  for  mt  =  4  and  m2  =  5. 

The  integers  within  each  modulus  are  representative  of  a  symmetrically  folded 
waveform  with  the  period  of  the  waveform  equal  to  the  modulus 

Psns  =  m  (6) 

where  m  is  the  modulus  (channel). 

The  dynamic  range  of  the  SNS  MSNS  <  Mosns  is  given  as 
(i)  For  PRP  moduli  m„  with  one  of  the  moduli  even,  then 
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where  j  ranges  from  1  to  N  -1  and  mj2,  m^,  ...  m,w  range  over  all  permutations  (2,  3, 
Nj.  The  product  of  j,  from  2  to  I  is  empty  and  its  value  is  1. 


(ii)  //  all  the  PRP  moduli  in  the  system  are  odd,  then 
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where  j  ranges  from  1  to  N  -  1  and  mu,  mi2,  mi3,  ...  miN  range  over  all  permutations  {  1,  2, 
3,...,NJ. 
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Figure  1.2:  The  Folding  Waveform  of  the  SNS  for  m,  -4  and  m2  =  5 


3.  Optimum  Symmetrical  Number  System 

The  optimum  symmetrical  number  system  (OSNS)  formulation  is  a  direct 
consequence  of  the  need  to  extract  the  maximum  amount  of  information  from  a 
symmetrically  folded  waveform  [7].  The  OSNS  scheme  is  also  composed  of  a  number  of 
PRP  moduli  m.  Integers  within  each  OSNS  modulus  is  defined  as 

xm=[0,l,...,m-l,m-l 1,0].  (9) 

Figure  1.3  shows  the  OSNS  folding  waveform  for  mj  =  4  and  m2  =  5.  The  integers 
within  each  SNS  modulus  are  representative  of  a  symmetrically  folded  waveform  with  the 
period  of  the  waveform  equal  to  the  twice  the  modulus  and 

P  osns  =  2m  (10) 

where  m  is  the  modulus.  Due  to  the  presence  of  ambiguities,  the  integers  within  (9)  do 
not  form  a  complete  system  of  length  2m  by  themselves.  The  ambiguities  that  arise 
within  the  modulus  are  resolved  by  considering  the  paired  values  from  all  moduli 
together.  By  recombining  the  N  moduli,  the  OSNS  is  rendered  a  complete  system  having 
a  one-to-one  correspondence  with  the  residue  number  system.  For  N  equal  to  the  number 
of  PRP  moduli,  the  dynamic  range  of  the  OSNS  is 


N 


MosNs=Tlmt  00 


i=i 


where  mt  is  the  i    modulus,  1  <i  <N  and  N  is  the  number  of  moduli.  The  dynamic  range 
is  also  the  position  of  the  first  repetitive  moduli  vector. 

4.         Application  for  Modular  Number  Systems  in  Signal  Processing 

A  modular  number  system  such  as  the  RNS  has  been  used  in  high  performance 
digital  signal  processing  [8].  Applications  include  a  high  rate  arithmetic  processor, 
digital  filters,  memory  compression,  and  image  processing.  Recently,  symmetrical 
number  systems  have  been  investigated  in  such  application  as  Analog-to-Digital 
Converters  (ADCs)  [7,  9-13],  Digital  Antennas  [14],  Undersampling  Receivers  [6,  15], 
and  Direction  Finding  Antennas  [14].  The  symmetrical  number  systems  offer  a 
considerable  amount  of  flexibility  in  implementation.     For  higher  resolution,  faster 


conversion  speed,  and  lower  power  dissipation  the  symmetrical  number  systems  play  an 
important  role  in  folding  ADCs  applications. 
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Figure  1.3:  The  Folding  Waveform  of  the  OSNS  for  mj  =  4  and  m2  =  5 


5.         Encoding  Errors  of  Modular  Number  Systems 

Each  of  the  above  modular  number  systems  for  which  the  integer  values  within 
each  modulus,  when  considered  together,  change  one  or  more  at  the  next  position  is 
shown  in  Table  1.1.  Consequently,  when  some  integers  at  a  position  to  change,  do 
change,  while  others  do  not,  there  will  be  a  large  error  (encoding  error).  Table  1.1  shows 
the  integer  values  within  the  modular  systems  for  moduli  mi  =  4  and  mj  =  5. 

When  the  input  is  1 ,  the  integer  values  within  the  modular  number  systems  are  1 
for  both  moduli.  By  increasing  the  input  from  1  to  2,  all  integer  values  within  the 
modular  number  systems  change  from  1  to  2.  If  one  of  the  integer  values  does  not 
change,  the  new  moduli  pairs,  for  example  (1,2)  or  (2,1)  correspond  to  new  input  values. 
For  the  RNS,  the  input  values  can  be  6  or  17.  For  the  SNS,  the  input  values  can  be  6,  13, 
14,  or  17.  And  for  the  OSNS,  the  input  values  can  be  17  or  18.  These  values  cause  large 
errors  (encoding  errors)  that  must  be  dealt  with.  These  large  encoding  errors  present  a 
number  of  difficulties  when  used  for  signal  processing.  The  goal  of  this  thesis  is  to 
develop  a  new  number  SNS  formulation  that  minimizes  the  encoding  errors.  To  examine 
the  performance  of  the  new  number  system  formulation  a  folding  ADC  architecture  is 
investigated. 
B.         PRINCIPAL  CONTRIBUTION 

A  new  symmetrical  number  system  with  applications  in  parallel  signal  processing 
is  investigated.  The  "Robust  Symmetrical  Number  System"  (RSNS)  is  a  modular  system 
in  which  the  integer  values  within  each  modulus,  when  considered  together,  change  one 
at  a  time  at  the  next  position  (Gray  code  properties).  Although  the  observed  dynamic 
range  of  the  RSNS  is  somewhat  less  than  the  optimum  symmetrical  number  system,  the 
Gray  code  properties  make  it  particularly  attractive  for  folding  analog-to-digital 
converters.  With  the  RSNS,  the  encoding  errors  (due  to  comparator  thresholds  not  being 
crossed  simultaneously)  are  eliminated,  as  is  the  need  for  the  corresponding  interpolation 
signal  processing  (reduced  complexity).  Computer  generated  data  is  used  to  help 
determine  the  properties  of  the  RSNS.  These  properties  include  the  largest  dynamic 
range  (number  of  distinct  consecutive  vectors),  and  the  position  of  the  largest  dynamic 


range  within  the  system.  The  position  of  the  maximum  unambiguous  dynamic  range  is 
also  quantified.  Least  squares  analysis  of  2  and  3  moduli  systems  is  used  to  formulate 
closed-form  expressions  for  the  dynamic  range.  To  compare  the  advantages  of  the  RSNS 
with  previously  published  results,  the  transfer  function  of  a  3-channel  RSNS  folding 
analog-to-digital  converter  architecture  (ray  =3,  m.2  -  4,  and  raj  =  5)  is  numerically 
evaluated  using  SPICE. 
C.         THESIS  OUTLINE 

This  thesis  is  organized  in  six  parts.  Chapter  II  provides  the  definition  of  the 
robust  symmetrical  number  system  (RSNS).  An  expression  for  the  fundamental  period  is 
given  and  proven.  In  this  chapter,  the  algorithm,  which  finds  the  dynamic  range,  point 
dynamic  range,  and  fundamental  period  is  discussed.  The  beginning  and  ending  of  the 
discrete  states  without  redundancy  for  any  given  number  of  channel  and  shift  values  of 
the  RSNS,  is  presented. 

Chapter  HI  investigates  the  2-channel  (moduli)  RSNS.  A  specific  number 
combination  is  considered.  The  RSNS  performance  using  different  moduli,  the  same 
distance  apart  is  evaluated  to  find  a  closed-form  expression  dynamic  range  in  the  general 
case. 

Chapter  IV  deals  with  the  3-channel  (moduli)  RSNS.  A  specific  number 
combination  is  also  considered  in  order  to  find  a  closed-form  expression  for  the  dynamic 
range. 

Chapter  V  is  based  on  the  simulation  of  a  3-channel  RSNS  ADC  in  SPICE.  The 
simulation  channel  and  the  shift  values  are  ray  =  3,  Sj  =  0,  mz  =  4,  S2  =  1,  and  raj  =  5,  sj 
=  2,  where  ra  shows  the  channel  values  and  s  shows  the  shift.  The  folding  circuit  and 
comparator  stages  are  designed  to  fold  and  quantize  the  incoming  signal  in  a  RSNS 
format. 

Chapter  VI  summarizes  the  results  and  conclusions  of  this  study  and  provides 
several  recommendations  for  follow-up  research. 
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Table  1.1:  Modular  Number  Systems  to  Decimal  Mapping  for  m;  =  4  and  m2  =  5 
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II.     ROBUST  SYMMETRICAL  NUMBER  SYSTEMS 


A.        DEFINING  THE  RSNS  WAVEFORM 

In  the  robust  symmetrical  number  system  (RSNS),  N  different  periodic 
symmetrical  waveforms  are  used  with,  pairwise  relatively  prime  (PRP)  integer,  lengths 
mi,  ni2,  ...,  mN.  The  RSNS  is  based  on  the  following  sequence 

Xmi  =  [  0,  1,  2 m  -  J,  m,  m-  J,...,  2,  1]  (12) 

where  xm]  is  the  row  vector  and  m  is  a  positive  integer  (m  >  0).  In  an  TV-channel  RSNS, 
where  N  >2,  the  basic  sequence  for  ith  channel,  with  modulus  m  is 

xm2=[0,  0,...,  0,  0,  1,  1,...,  1,  1,...,  m,  m,...,  m,  m  ....  /,  1,...,  1,  I].  (13) 

In  this  sequence  each  value  in  the  xmj  row  vector  is  put  TV  times  in  succession.  This 
sequence  is  repeated  in  both  directions,  forming  a  periodic  sequence  with  the  period 

PRSNS  =  2mN.  (14) 

An  /V-channel  RSNS  is  formed  of  vectors  by  picking  TV  moduli  (channels)  Wj,  and  N  shift 
values  su  1  <i  <N,  where  {si,  S2,...,  s^}  forms  a  complete  set  of  residues  modulo  N.  The 
Vh  vector 


Vh  = 


■i/i 


a 


m 


(15) 


■th 


where  an,  is  the  (h  -  Sj)  term  in  the  periodic  sequence  for  w„  1  <i  <N.  Since  the  /  row  is 
periodic  with  period  2Nm„  1  <i  <N,  it  is  clear  that  if  the  length  that  is  simultaneously  a 
multiple  of  2Nmi,  2Nm2,...,  and  2NmN  is  taken,  then  the  N-channel  RSNS  vector  will 
repeat  itself.  Since  the  fundamental  period  for  channel  i  is  2Nm„  J  <i  <N,  it  follows  that 
a  period  for  the  RSNS  vectors  must  be  a  multiple  of  each  of  2Nm,.  Therefore  the 
fundamental  period  for  the  RSNS  is 

PFrsns  =  [2m jN,  2m2N, ....  2mNN]  ( 1 6) 
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where  [nj,  n2,  ...,  /W  is  the  least  common  multiple  of  nj,  n2,  ...,  n^.  From  number  theory 
[2-5], 

PFrsns  =  2N[mi,  m2,...,  mN].  ( 1 7) 

Considering  a  single  channel  the  discrete  states  of  the  robust  symmetrical  number 
system  are  given  as 

n  —  5. 


g(mi,si,n)  =  - 


2Nm:  +  N-n  +  si-\ 
N 


for  1  +  st <  n  <  Nmi  +  st 


for  Nmi  +  \  +  sj  <n<  2Nmj  +  st 


(18) 


where  m,  is  the  channel  modulus  and  Sj  is  the  corresponding  shift  for 

Si=0,l,2,...,N-l(modN)  (19) 

and  N  >2  is  the  number  of  channels  in  the  system.  The  values  {sj,  S2,...,  W  must  form  a 
complete  residue  system  modulo  N.  Here  n  represents  the  normalized  input  voltage. 
Because  of  the  relative  property  of  the  shifts,  one  of  the  shift  values  will  be  set  equal  to  0. 

Upon  inspection  of  any  2-channel  and  3-channel  RSNS,  it  is  seen  that  if  Vh  (the 
vector  at  position  h)  equals  Vh+k,  then  k  is  a  multiple  of  4  and  6,  respectively.  These  are 
particular  examples  of  a  general  rule: 

Rule  2.1:  Suppose  that  Vh,  Vh+k  cire  the  vectors  at  position  h  and  h  +  k,  respectively,  in  an 
N -channel  RSNS.  IfVh  =  Vh+k,  then  k  is  the  multiple  of2N. 

The  proof  is  given  below,  but  first  recall  that  in  an  TV-channel  RSNS  each  term  in  a 
given  channel  occurs  consecutively  N  times.  Consider  the  TV  consecutive,  equal  terms 
form  from  left  to  right  as  lying  in  the  first,  second,. . .,  Nth  position.  Since  2N  is  a  multiple 
of  N,  this  rule  shows  that  at  any  redundancy  the  redundant  term  lies  in  the  same  relative 
position  as  the  original  term.  That  is,  if  Vh  =  Vh+k, 

a, 


Vh  = 


lNh 


>    Vh+k 


ll.h  +  k 


lN.h+k 


(20) 
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and  aAj,  occurs  in  the  i    position  in  row  A.   Then,  aA,h  =  aA.h+k  and  aAj,+k  also  occurs  in 
the  i'h  position  in  row  A,  1  <A  <N. 

The  proof  of  the  Rule  2.1: 

In  the  fundamental  building  block  of  the  RSNS,  the  sequence  by  (12),  the  value  at 
j  (0  <j  <  m)  equals  to  the  value  at  (j  +  k)  when  j  +  k  =  m  +  (m  -j),  i.e.  when  k  =  2m  -  2j, 
an  even  number.  This  fundamental  building  block  has  length  2m,  and  it  is  extended 
periodically.  The  values  0,  and  m  repeat  every  2m  terms,  and  no  more.  Other  values 
repeat  twice  as  often,  but  always  an  even  number  of  terms  apart.  When  the  sequence  is 
expanded  by  the  TV-fold  repetition  of  each  term,  then  all  terms  that  repeat  in  the  same 
relative  position  are  a  multiple  of  2N  apart. 

Suppose  that  ah  =  ah+k,  k  >  0,  and  that  ah  is  in  the  first  (left)  position.  Then  ah+k 
will  be  just  i  terms  to  the  right  of  the  term  in  the  first  position,  where  0  <  i  <  N.  By  the 
conclusion  of  the  last  paragraph  k  =  2Nj  +  i  for  some  integers  i,  0  <i  <  N,  and  j  >  0  (  and 
not  both  0).  Similarly,  if  ah  =  ah+k,  where  ah  is  in  the  N*h  (right)  position,  then  it  follows 
that  k  =  2Nj  -  i ,  for  some  i  ,0  <i  <  N,  and  j  >  0. 

Now  return  to  the  full  A^-channel  problem.  Suppose  Vh  =  Vh+k  with  k  >  0.  By  the 
definition  of  the  RSNS  one  of  the  terms  in  (15)  is  in  the  first  position,  say  aAh,  and  one  is 
in  the  TV^  position,  say  a^.h-  Since  Vh+k  equals  Vh,  aA,h  =  tf/u+*  and  as.h  =  QB.h+k-  The  first 
of  these  equations  implies  that  k  =  2Nj  +  i  for  some  i,0  <  i  <  N  andy  >  0,  and  the  second 
equation  implies  that  k  =  2Nj  -  i ,  for  some  i ,  0  <  i  <  N,  and  j  >  0.  Putting  these 
together:  Since  2Nj  +  i  and  2Nj  -  i  are  both  equal  to  k,  2Nj  +  i  =  2Nj  -  i .  Thus  i  +  i  = 
2N(j  -j).  However,  0  <i  +  i  <  2N,  so  /  +  /  =  0  and  j  =  j  .  Since  /  and  i  are  nonnegative 
both  must  be  zero.  Therefore  k  is  a  multiple  of  2N,  as  claimed. 

The  discrete  states  of  the  RSNS  are  indexed  starting  from  the  vector,  which  has  all 
zero  values.  The  index  value  starts  from  0.  These  index  values  are  put  in  a  row  matrix  h, 
which  lies  in  {..., -2, -1,  0,  1,  2,  3,...}  and  these  values  are  called  the  point  index.  For 
2-channel  RSNS,  all  point  indexes  determine  2  non-negative  integer  numbers,  which  are 
the  discrete  states  of  the  RSNS.  Similarly,  //-channel  RSNS  point  indexes  determine  N 
discrete  states  of  the  RSNS. 
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Figure  2.1  shows  the  2-channel  RSNS  folding  waveforms  and  the  discrete  states 
within  the  RSNS.  The  thresholds  shown  on  the  vertical  axis  represent  the  integer  value 
within  each  RSNS  modulus.  The  integer  values  occur  2  times  in  succession.  In  this 
example  the  channel  and  corresponding  shift  values  are  my  =  4,  si  =  0  and  m.2  =  5,  S2  =  1. 

Figure  2.2  shows  3-channel  RSNS  folding  waveforms  and  the  discrete  states 
within  the  RSNS.  The  thresholds  shown  on  the  vertical  axis  represent  the  integer  value 
within  each  RSNS  modulus.  The  integer  values  occur  3  times  in  succession.  In  this 
example  the  channel  and  corresponding  shift  values  are  mj  =  3,  si  =  0,  m.2  -  4,  S2  =  1 
and  raj  =  5,  S3  =  2. 

To  demonstrate  the  performance  of  the  2  and  3-channel  RSNS  in  detail,  the 
RSNS  computer  algorithm  is  executed  and  Tables  2.1  and  2.2  show  the  distinct  vector 
positions.  Each  row  is  considered  separately.  The  first  and  second  column  of  a  given 
row  give  the  beginning  and  ending  point  indexes  containing  consecutive  sequences  of 
distinct  vectors  of  iV  states.  Furthermore,  these  positions  can  not  be  extended  without 
introducing  an  ambiguity.  The  total  length  of  the  consecutive  sequence  of  distinct 
vectors,  starting  at  the  beginning  points  (the  point  dynamic  range)  is  shown  in  the  third 
column.  Note  the  maximal  number  of  distinct  vectors  in  Table  2.1  is  21  and  in  Table  2.2 
is  43. 
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Beginning  Point 
Index 

Ending  Point 
Index 

Total 
Length 

1 

18 

18 

4 

19 

16 

5 

25 

21 

7 

26 

20 

8 

28 

21 

14 

29 

16 

15 

33 

19 

31 

42 

12 

Table  2.1:  The  Position  of  Distinct  Vectors  for  the  2-Channel  Case  RSNS 

(mi  =  4,  sj  =  0,  and  rri2  =  5,  S2  =  1) 


Beginning  Point 
Index 

Ending  Point 
Index 

Total 
Length 

2 

29 

28 

7 

30 

24 

8 

41 

34 

13 

51 

39 

23 

52 

30 

24 

57 

34 

41 

65 

25 

61 

103 

43 

Table  2.2:  The  Position  of  Distinct  Vectors  for  the  3-Channel  Case  RSNS 

(mi  =  3,  Sj  —  0,  rri2  =  4,  S2  =  1,  and  Wj  =  5,  sj  -  2) 
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B.       THE  DYNAMIC  RANGE  OF  AN  RSNS 

The  dynamic  range  of  the  RSNS  is  the  maximum  number  of  distinct  vectors 

without  a  redundancy.  The  dynamic  range  of  the  RSNS  is  shown  with  the  symbol  AtfMMS  . 
The  point  dynamic  range  of  the  RSNS  is  the  maximum  number  of  consecutive  vectors 
without  a  redundancy  starting  at  a  given  vector  (point  index).  The  point  dynamic  range  of 

the  RSNS  is  shown  with  the  symbol  MPSSNS  ■ 

The  performance  of  a  2-channel  RSNS  is  shown  in  Table  2.3  for  mj  =  4,  s}  =  0 
and  rri2  =  5,  s2  =  1.  The  selection  of  the  shifts  and  their  permutations  among  the  2-moduli 
has  no  effect  on  the  dynamic  range.  The  point  index  (beginning-ending)  of  the  set  of 
vectors  that  gives  the  dynamic  range,  however,  is  different. 


Moduli 
(mi) 

Shift 

Dynamic 
Range 

V  *™ RSNS  ' 

Fundamental 
Period 

(Prsns) 

Beginning  -  Ending  Point 

Index  of 

Dynamic  Range 

4 
5 

0 

i 

21 

80 

5-25 

4 
5 

i 

0 

21 

80 

14-34 

Table  2.3:  The  Performance  with  the  Permutation  of  Shifts  for  the  2-Channel  Case 

RSNS 

The  performance  of  a  3-channel  RSNS  is  shown  in  Table  2.4  for  m;  =  3,  sj  =  0, 
n%2  -  4,  s2  =  1  and  m3  =  5,  S3  =  2.  The  selection  of  the  shifts  and  their  permutations 
among  the  3-moduli  also  has  no  effect  on  the  dynamic  range.  The  point  index 
(beginning-ending)  of  the  set  of  vectors  that  determines  the  dynamic  range  is  different. 
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Moduli 

Shift 

Dynamic 

Fundamental 

Beginning  -  Ending  Point 

(md 

(Si) 

Range 

Period 

Index  of 

V  M  RSNS  ' 

(Prsns) 

Dynamic  Range 

3 

0 

4 

1 

43 

360 

61-103 

5 

2 

3 

0 

4 

2 

43 

360 

142-184 

5 

1 

3 

1 

4 

0 

43 

360 

8-50 

5 

2 

3 

1 

4 

2 

43 

360 

134-176 

5 

0 

3 

2 

4 

0 

43 

360 

180-222 

5 

1 

3 

2 

4 

1 

43 

360 

81-123 

5 

0 

Table  2.4:  The  Performance  with  the  Permutation  of  Shifts  for  the  3-Channel  Case 

RSNS 
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C       A  WEAKNESS  OF  THE  RSNS  WITH  RESPECT  TO  THE  INCREASE  OF 
THE  NUMBER  OF  CHANNELS 

The  integers  fall  into  two  classes,  even  and  odd  parities.  The  integer  values  for  any 
number  of  RSNS  channels  examined  by  a  parity  point  of  view.  For  the  2-channel  ideal 
case,  their  parities  distribute  22  =  4  different  ways.  These  are  even-even,  even-odd,  odd- 
even,  and  odd-odd.  The  possible  parity  sets  for  the  2-channel  RSNS  are  shown  in  Table 
2.5  and  these  parity  sets  repeat  themselves.  In  the  2-channel  RSNS,  the  efficiency  of  the 
parity  sets  is  same  as  the  ideal  one. 


Index 

mi 

m.2 

1 

0 

0 

2 

0 

1 

3 

1 

1 

4 

1 

2 

5 

2 

2 

6 

2 

3 

7 

3 

3 

8 

3 

4 

9 

4 

4 

10 

4 

5 

11 

12 

13 

Even-Odd 

Odd-Odd 

Odd-Even 

Even-Even 

Even-Odd 

-<J 

Odd-Odd 

Odd-Even 

Even-Even 

Table  2.5:  Demonstration  of  the  parity  combinations  for  N  =  2  with  si  =  0  and  s2  =  1 


For  the  3-channel  ideal  case,  their  parities  distribute  23  =  8  different  ways.  These 
are  even-cven-even,  even-even-odd,  even-odd-odd,  even-odd-even,  odd-even-even,  odd- 
even-odd,  odd-odd-even,  and  odd-odd-odd.   There  are  only  6  possible  parity  sets  for  the 
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3-channel  RSNS  (restricted  parity  combinations)  and  they  repeat  themselves  as  shown  in 
Table  2.6.  The  parity  sets  even-odd-even  and  odd-even-odd  will  never  show  up  in  this  3- 
channel  RSNS.  In  the  3-  channel  RSNS,  the  efficiency  of  the  parity  sets  is  25%  less  than 
the  ideal  one. 


Index 

m, 

m2 

m3 

1 

0 

0 

0 

2 

0 

0 

1 

3 

0 

1 

1 

4 

1 

1 

1 

5 

1 

1 

2 

6 

1 

2 

2 

7 

2 

2 

2 

8 

2 

2 

3 

9 

2 

3 

3 

10 

3 

3 

3 

11 

3 

3 

4 

12 

3 

4 

4 

13 

14 

15 

Even-Even-Even 

Even-Even-Odd 

Even-Odd-Odd 

Odd-Odd-Odd 

■«< 

Odd-Odd-Even 

Odd-Even-Even 

Even-Even-Even 

Even-Even-Odd 

■^ 

Even-Odd-Odd 

Odd-Odd-Odd 

<4 

Odd-Odd-Even 

Odd-Even-Even 

Table  2.6:  Demonstration  of  restricted  parity  combinations  for  N  =  3  with 

sj  =  0,  S2  =  1  and  S?  =  2 


In  the  TV-channel  RSNS  each  term  in  the  vector  Vh+N  by  (20)  differs  by  exactly  one 
from  corresponding  term  in  V/,  by  (15).  Therefore,  each  term  in  Vh+N  has  opposite  parity 
from  the  corresponding  term  in  V/,.  It  follows  quickly  that  each  term  in  Vh+2N  has  exactly 
the  same  parity  as  the  corresponding  term  in  Vh.   From  this  result,  it  can  be  said  that  the 
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parities  of  the  terms  in  any  V/,  will  be  exactly  the  same  as  the  parities  of  the  corresponding 
terms  in  one  of  Vj,  V2,  ...,  V^.  This  conclusion  ends  the  proof  that  there  will  be  just  2N 
different  parity  sets  represented,  out  of  a  possible  2N  for  ^/-channel  RSNS. 
D.       RSNS  COMPUTER  ALGORITHM 

In  this  study,  one  of  the  important  goals  is  to  be  able  to  find  a  closed-form 
expression  for  the  dynamic  range  and  fundamental  period  of  any  number  of  RSNS 
channels.  The  fundamental  period  of  any  number  of  RSNS  is  given  in  Section  A.  The 
RSNS  computer  algorithm  (rsns.m)  is  constructed  in  order  to  examine  the  performance 
of  the  RSNS  and  find  the  dynamic  range  for  any  N  >  2  RSNS  systems. 

In  rsns.m,  the  channel  and  shift  values,  which  are  going  to  be  used  must  be  entered 
and  the  discrete  states  of  the  RSNS  are  calculated  by  using  these  values  along  with  (18). 
The  calculated  discrete  states  of  the  RSNS  are  put  into  a  matrix.  This  matrix  is  extended 
two  times  the  fundamental  period  of  the  RSNS  channels. 

The  next  step  in  the  algorithm  is  to  find  the  point  index  of  the  discrete  states  within 
the  extended  matrix.  After  finding  the  point  index,  the  data  is  put  into  a  new  matrix.  The 
point  index  of  each  discrete  state  and  its  redundancy  is  assigned  as  the  beginning  and 
ending  point  index,  respectively. 

The  values  of  the  new  matrix  are  sorted  and  rearranged  such,  that  the  beginning  and 
ending  point  index  are  monotonically  increasing.  Then,  the  length  at  each  point  index  is 
found  (Ending  point  index  -  Beginning  point  index  +  1).  The  maximum  value  of  these 
lengths  is  the  dynamic  range  of  the  RSNS. 

The  algorithm  rsns.m  also  displays  the  dynamic  range,  fundamental  period,  entered 
channel  and  shift  values,  the  beginning-ending  point  index  values  of  the  discrete  states 
without  redundancy  and  the  point  dynamic  range  of  the  RSNS.  The  program  structure 
and  listing  of  the  rsns.m  is  shown  in  Appendix.  Tables  2.1  and  2.2  are  the  example  runs 
of  the  rsns.m  algorithm. 
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III.  2-CHANNEL  ROBUST  SYMMETRICAL  NUMBER  SYSTEM 

A.        SPECIAL  CASE:  MODULI  SPACING  1  OR  2  APART 

In  this  chapter,  the  2-channel  RSNS  is  addressed.  The  2  moduli  are  RP  numbers 
and  the  different  pairs  are  chosen  to  have  the  same  distance  between  each  other.  The 
algorithm  rsns.m  is  used  to  find  the  dynamic  range.  When  the  moduli  have  a  difference 
(A)  of  1  or  2,  the  following  computer  results  have  been  observed  as  shown  in  Table  3.1 
and  3.2. 


nt] 

m2 

Dynamic  Range 

1st  Difference 

3 

4 

15 

6 
6 
6 
6 

4 

5 

21 

5 

6 

27 

6 

7 

33 

Mksns 

=  3mi  +  3m2  -  6 

Table  3.1:  The  dynamic  range  of  mj  and  m2  =  mi  +  J  where  mi  >  3  (A  =  1) 


m2 

mi 

Dynamic  Range 

1st  Difference 

1 

5 

29 

12 
12 
12 
12 

9 

7 

41 

11 

9 

53 

13 

11 

65 

Mrsns  -  3m!  +  3m2  -  7 

Table  3.2:  The  dynamic  range  of  my  and  m2  =  mi  +  2  where  my  >5  (A  =  2) 


These  relationships  show  a  uniform  increase  of  3  with  each  unit  increase  of  value 
in  my  and  m2.  The  formulas 
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M 


RSNS 


M 


RSNS 


=  3mj  +  3rri2  -  6 
=  3mj  +  3rri2  -  7 


(21) 
(22) 


will  work  for  the  first  and  second  sets  of  cases  which  are  shown  in  Table  3.1  and  3.2, 
respectively.  These  cases  can  be  expressed  in  terms  of  the  lower  value,  my  and  a 
difference,  A-m2-mi.  In  the  first  set  A  =  1,  and  in  the  second  set  A  =  2.  Thus,  (21)  and 
(22)  can  be  written  as 


M 


RSNS 


M 


RSNS 


=  3mi  +  3(mj  +  A)  -  6  =  6mj  -  3 
=  3mj  +  3(mj  +  A)  -  7  =  6m\  -  1. 


These  can  be  combined  as 


M 


RSNS 


=  6mi  +  2A-  5 


and  now  A  can  be  replaced  with  ni2  -  mi,  yielding 


M 


RSNS 


=  6mi  +  2(m2-mi)  -  5  =  4mi  +  2m2  -  5 


(23) 
(24) 

(25) 

(26) 


and  this  formula  works  for  both  A  =  1  and  A  =  2.   This  formula  is  checked  against  the 
computer  generated  results.  For  example,  when  mi=J00  and  m2=10J,  the  dynamic  range 
is  597  which  fits  (26). 
B.        SPECIAL  CASE:  MODULI  SPACING  3  OR  MORE  APART 

The  following  results  have  been  observed,  when  the  channel  numbers  are  3,  4  and 
5  distance  apart  as  shown  in  Table  3.3,  3.4  and  3.5. 


mi 

m2 

Dynamic  Range 

1st  Difference 

5 

8 

34 

18 

18 
18 
18 

8 

11 

52 

11 

14 

70 

14 

17 

88 

Mrsns  =  3m>  +  3m2  -  5 

Table  3.3:  The  dynamic  range  of  my  and  ni2  =  mj  +  3  where  my  >  5  (A  =  3) 
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mj 

m2 

Dynamic  Range 

1st  Difference 

1 

11 

48 

24 
24 
24 
24 

11 

15 

72 

15 

19 

96 

19 

23 

120 

M 

lrx  RSNS 

=  3mi  +  3m2  -  6 

Table  3.4:  The  dynamic  range  of  mj  and  m2  =  ray  +  4  where  mj  >7  (A-  4) 


mi 

m2 

Dynamic  Range 

1st  Difference 

1 

12 

50 

30 
30 
30 
30 

12 

17 

80 

17 

22 

110 

22 

27 

140 

™  RSNS 

=  3mj  +  3m2  -  7 

Table  3.5:  The  dynamic  range  of  my  and  m2  =  mj  +  5  where  m;  >  7  (A  =  5) 

These  relationships  show  a  uniform  increase  of  3  with  each  unit  increase  of 
value  in  mt  and  m2.    The  formulas 

MRSNS  =  3mi  +  3m2  -  5  (27) 

M rsns  ~  3mi  +  3m2  -  6  (28) 

Mrsns    ~  3mj  +  3m2  -  7  (29) 

will  work  for  the  first,  second  and   third  sets  of  cases  which  are  shown  in  Table  3.3,  3.4 

and  3.5.  respectively.. 

As  in  the  earlier  example  these  cases  can  be  expressed  in  terms  of  the  lower  value, 

mi  and  a  difference,  A  =  m2  -  mi.  In  the  first  set  A  =  3,  in  the  second  set  A  =  4,  and  in  the 

third  set  A  =  5.  Thus,  (27),  (28)  and  (29)  can  be  written  as 
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MRSNS  =  Smi  +  3(mi  +  A)  -  5  =  6m1  +  4  (30) 

MRSNS  =  3m i  +  3(mi  +  A)  -  6  =  6mj  +  6  (31) 

Mrsns  =  3mi  +  3(mi  +  A)  -  7  =  6mi  +  8.  (32) 

These  can  be  combined  as 

MgsNs  =6rm  +  2A-2  (33) 

and  now  A  can  be  replaced  with  m.2  -  mi,  yielding 

MRSNS  =  6mj  +  2(m2-mi)  -  2  =  4mi  +  2m2  -  2.  (34) 

and  this  formula  works  for  all  three  values  of  A,  A  =  3,  A  =  4,  and  A  =  5.  This  formula  is 
checked  against  the  computer  generated  results.  For  example,  when  mi=J0J  and 
m2=105,  the  dynamic  range  is  612  which  fits  (34).  It  is  conjectured  that  (34)  is  valid  for 
any  A  >  3.  Table  3.6  shows  a  summary  of  the  specific  2-channel  dynamic  range  formulas. 


Spacing  between  moduli  (A  ) 

Dynamic  Range  Formula 

A=l  or  A  =  2 

MRSNS    =  4m  i  +  2mj  -  5 

A>3 

MRSNS  =  4m  i  +  2m2  -  2 

Table  3.6:  The  Summary  of  the  Dynamic  Range  Formulas  for  specific  2-channel 

RSNS 

C.         GRAPHICAL  ANALYSIS 

The  algorithm  rsns.m  also  finds  and  displays  the  point  dynamic  range  of  the  2- 
channel  RSNS  as  explained  in  Chapter  H  The  point  dynamic  range  versus  the  point 
index  for  the  channel  and  shift  values  m}  =  4,  S]  =  0  and  m2  =  5,  s2  =  1  are  shown  in 
Figures  3.1  and  3.2,  respectively.   The  waveform  with  the  points  A,  B,  C,  D  and  E,  that 

p 

are  shown  in  Figures  3.1  and  3.2,  construct  a  special  waveform  that  repeats  itself     RSNS  . 

This  special  waveform,  the  point  dynamic  range  (PDR)  waveform,  gives  the  dynamic 
range  starting  at  each  point  with  increasing  point  index.   The  fundamental  period  of  this 
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particular  2-  channel  RSNS  example  is  PFrsns  =  80  by  (16).  So,  this  PDR  waveform 
repeats  itself  every  40  points. 

The  maximum  value  of  the  point  dynamic  range  is  the  dynamic  range  of  the  2- 

channel  RSNS.  In  this  example,  MRSNS  =21.  The  dynamic  range  point  indexes  are  5,  8, 
45,  48,  ...  for  si  =  0,  S2  =  1  and  14,  17,  54,  57,  ...  for  si  =  1,  S2  =  0  which  can  be  seen 
from  Figure  3.1  and  3.2.  These  values  are  the  beginning  point  index  of  the  dynamic 
range.  The  dynamic  range  point  index  values  are  also  half  of  the  fundamental  period 
apart  from  each  other. 

The  PDR  waveform  of  the  point  dynamic  range  of  the  2-channel  RSNS  will  be  the 
same  for  different  shift  values  as  shown  in  Figure  3.1  and  3.2.  The  only  difference  is  the 
point  index  (beginning  value)  of  the  dynamic  range. 
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40  50  60 

The  Point  Index 


100 


Figure  3.1:  The  Point  Dynamic  Range  Waveform  for  nij  —4,  sj  =  0,and  m2  =  5, 

s2  =  l. 
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Pomi 

Point  Index 

Mr.KSK 

A 

9 

19 

B 

14 

21 

C 

17 

21 

D 

24 

19 

E 

40 

12 

40  50  60 

The  Point  Index 


100 


3.2:  The  Point  Dynamic  Range  Waveform  for  m\  -  4,  s\-  7, and  m.2  =  5, 

s2  =  0. 
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IV.      3-CHANNEL  ROBUST  SYMMETRICAL  NUMBER  SYSTEM 

A.         SPECIAL  CASE:  MODULI  SPACING  1  APART 

In  this  chapter,  the  3-channel  RSNS  is  addressed.  An  equal  distance  separates  the 
moduli.  When  the  PRP  channel  numbers  are  1  apart  and  the  set  of  the  3-channel  RSNS 
are  increasing  uniformly,  the  following  results  have  been  observed  as  shown  in  Table  4. 1 . 


mi 

rri2 

rri3 

Dynamic  Range 

1st  Difference 

2nd  Difference 

3 

4 

5 

43 

39 

51 
63 

12     • 
12 

5 

6 

7 

82 

7 

8 

9 

133 

9 

10 

11 

196 

3    ,     15 

MRSNS   =   2mi'+Ymi+7 

Table  4.1:  The  Dynamic  Range  of  my,  m.2  =  my  +  7,  and  m3  =  my  +  2  where  my  >3 

The  first  difference  (e.g.,  82  -  43  =  39)  and  second  difference  (e.g.,  51  -  39  =  12) 
of  the  dynamic  range  of  the  3-channel  RSNS  for  mj,  m2  =  mi  +  1,  and  mj  =  my  +  2  where 
mi  >  3  are  shown  in  Table  4. 1 .  The  relationship  between  the  dynamic  range  differences 
can  be  expressed  as  a  2nd  degree  polynomial 

p(n)  =  an2  +  bn  +  c  (35) 

where  n  is  the  lowest  channel  value  and  a,  b,  and  c  are  a  constant  real  numbers.  Each  set 
of  3-moduli  RSNS  values  shows  a  uniform  increase  of  2.  Therefore,  if  the  first  set  of 
channel  values  corresponds  to  (35),  then  the  second  set  of  channel  values  corresponds  to 

p(n+2)  =  a(n  +  2)2  +  b(n  +  2)  +  c  (36) 

=  a(n2  +  4n  +  4)  +  bn  +  2b  +  c 
=  an2  +  4an  +  4a  +  bn  +  2b  +  c 
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=  p(n)  +  4an  +  4a  +  2b. 
The  first  difference  of  the  sets  ((36)  -  (35))  is 

q(n)  =  p(n  +  2)  -  p(n)  =  4an  +  4a  +  2b.  (37) 

If  the  first  difference  of  the  first  two  sets  is  (37),  the  first  set  difference  of  the  second  two 
sets  can  be  written  as 

q(n  +  2)  =  4a(n  +  2)  +  4a  +  2b  (38) 

=  4an  +  8a  +  4a  +  2b 
=  q(n)  +  8a. 
The  second  difference  of  the  sets  are  unique  and 

q(n+2)  -  q(n)  =  8a.  (39) 

The  constant  "a"  can  be  found  by  using  (39) 

q(n+2)  -  q(n)  =  8a  =12  (40) 

\2_  3 

a~  8  ~2 

The  constant  "b"  can  be  found  by  using  (37)  for  n  =  3  (  mi  =  3) 

f3\         f 


p(n  +  2)  -  p(n)  -  39  =  4 


-  3  +  4 


-\  +  2b 
2 


36  +  12 
39  =  — - —  +  2b  (41) 

2 

15 

h  =  Y 

The  constant  "c"  can  be  found  by  using  (35)  for  n  =  3  (mj  =  3) 

3    ,     15  3   ,     15 

p(3)  =  43  =  -n2+—n  +  c  =  -32+—3  +  c 
y  2  2  2  2  (42) 

c  =  l 
Finally,  the  closed-form  expression  for  the  dynamic  range  can  be  conjectured  as 

3    ,      15 
Mrsns  =  2mi'+Ymi+7  (43) 

where  m,  >  3.  This  formula  is  checked  against  the  computer  results.   For  example,  when 
mi  =  23,  m2  =  24,  and  mj  =  25,  the  dynamic  range  is  973  which  fits  (43). 
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B.         SPECIAL  CASE:  MODULI  SPACING  2-  APART 

When  the  PRP  channel  numbers  are  2  apart,  the  following  results  have  been 
observed  as  shown  in  Table  4.2. 


ray 

m.2 

m3 

Dynamic  Range 

1st  Difference 

2nd  Difference 

5 

1 

9 

104 

126 

174 
222 

48 
48 

9 

11 

13 

230 

13 

15 

17 

404 

17 

19 

21 

626 

3    ,     21 
Mrsns  =  2^  +y^+14 

Table  4.2:  The  Dynamic  Range  of  ray,  m2  =  mj  +  2,  and  raj  =  ray  +  4  where  mi  >  5 


The  first  difference  (e.g.,  230  -  104  =  126)  and  second  difference  (e.g.,  174  -  126 
=  48)  of  the  dynamic  range  of  the  3-channel  RSNS  for  ra/,  ra2  =  ray  +  2,  and  m3  =  mi  +  4 
where  ra/  >  5  are  shown  in  Table  4.2.  The  relationship  between  the  dynamic  range 
differences  can  also  be  expressed  as  a  2n  degree  polynomial 

p(n)  =  an2  +  bn  +  c  (44) 

where  n  is  the  lower  channel  value  and  a,  b,  and  c  are  constant  real  numbers.  Each  set  of 
the  3-moduli  RSNS  values  shows  a  uniform  increase  of  4.  Therefore,  if  the  first  set  of 
channel  values  corresponds  to  (44),  then  the  second  set  of  channel  values  corresponds  to 


p(n+4)  =  a(n  +  4)   +  b(n  +  4)  +  c 

=  a(n    +  8n+  16)  +  bn  +  4b  +  c 
=  an2  +  8an  +  16a  +  bn  +  4b  +  c 
=  p(n)  +  San  +  16a  +  4b. 
The  first  difference  of  the  sets  ((45)  -  (44))  is 

q(n)  =  p(n  +  2)  -  p(n)  =  8an  +  16a  +  4b. 
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(45) 


(46) 


If  the  first  difference  of  the  first  two  sets  is  (46),  the  first  set  difference  of  the  second  two 

sets  can  be  written  as 

q(n  +4)  =  8a(n  +  4)+  16a  +  4b  (47) 

=  8an  +  32a  +16a  +  4b 

=  q(n)  +  32a. 

The  second  difference  of  the  sets  are  unique  and 

q(n+2)  -  q(n)  =  32a.  (48) 

The  constant  "a"  can  be  found  by  using  (48) 

q(n+4)  -  q(n)  =  32a  =48  (49) 

48      3 
a=V  =  2 

The  constant  "b"  can  be  found  by  using  (46)  for  n  =  5  (mi  =  5) 


p(n  +  4)-/7(n)  =  126  =  8(-j5+l 


-  \  +  4b 


120  +  48 
126  = +  4/7  (50) 

2 

21 
b  =  — 

2 

The  constant  "c"  can  be  found  by  using  (44)  for  n  =  5  (ntj  =  5) 

3         21  3         21 

p(5)  =  104  =  -n2  +—  n  +  c  =  -52  +—  5  +  c 

2  2  2  2  (51) 

c  =  14 
Finally,  the  closed-form  expression  for  the  dynamic  range  can  be  conjectured  as 

3  21 

MRSNS  =  2ml+Ym]  +  l4  (52) 

where  m/  >  3.  This  formula  is  checked  against  the  computer  results.  For  example,  when 
mi  =  25,  m2  =  27,  and  raj  =  29,  the  dynamic  range  is  1214  which  fits  (52). 
C.         GRAPHICAL  ANALYSIS 

The  algorithm  (rsns.m)  also  finds  and  displays  the  point  dynamic  range  of  the  3- 
channel  RSNS  as  explained  in  Chapter  EL  The  point  dynamic  range  versus  the  point 
index  for  the  channel  values  mi  =  3,  m2  =  4,  and  =  5  are  shown  in  Figures  4.1  through 
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4.6.  The  waveform  with  the  points  A,  B,  C,  D,  E,  F,  G,  H,  I,  J,  K,  L,  and  M,  that  are 
shown  in  Figures  4. 1  through  4.6,  construct  a  special  waveform  that  repeats  itself  every 

p 

.     This  special  waveform,  the  point  dynamic  range  (PDR)  waveform,  gives  the 

dynamic  range  starting  at  each  point.  The  fundamental  period  of  this  particular  3-channel 
RSNS  example  is  360  by  (16).  So,  this  spike  waveform  repeats  itself  every  180  points. 

The  maximum  value  of  the  point  dynamic  range  is  the  dynamic  range  of  the  3- 
channel  RSNS.  In  this  example,  the  dynamic  range  is  43.  The  first  point  index  that 
yields  the  dynamic  range  is  61,  142,  8,  134,  180,  and  81  for  Figures  4.1  through  4.6, 
respectively.  These  values  are  the  beginning  point  index  of  the  dynamic  range.  The 
dynamic  range  point  index  values  are  also  half  of  the  fundamental  period  apart  from  each 
other. 

The  PDR  waveform  of  the  point  dynamic  range  of  the  3-channel  RSNS  will  be  the 
same  for  different  shift  values  as  shown  in  Figures  4. 1  through  4.6.  The  beginning  value 
of  the  dynamic  range  is  different  from  each  other.  When  one  of  the  shift  values  is  fixed, 
the  other  two  shift  value  combinations  of  the  PDR  waveforms  are  opposite  of  each  other. 
The  PDR  waveform  of  the  Figures  4.1,  4.3,  and  4.5  are  the  opposite  of  Figures  4.2,  4.4, 
and  4.6,  respectively. 


35 


40 


35 


30 


£20 

CD 


15 


10 


K 


G 


D 


H 


Poinl 

Poinl Index 

**r«j« 

A 

61 

43 

B 

81 

31 

C 

107 

24 

D 

IIS 

34 

E 

126 

34 

F 

132 

33 

G 

143 

35 

H 

161 

24 

I 

180 

30 

J 

188 

34 

K 

193 

39 

L 

204 

34 

M 

221 

25 

50  100  150  200  250 

The  Point  Index 


300 


350 


400 


Figure  4.1:  The  Point  Dynamic  Range  Waveform  for  mj  =  3,  sj  =  0,  m.2  =  4,  S2  -  1, 

m3  =  5,  s3  =  2. 
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Figure  4.2:  The  Point  Dynamic  Range  Waveform  for  mj  =  3,  sj  =  0,  m2  =  4,  s2  =  2, 

m.3  =  5,  S3  -  I . 
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Figure  4.3:  The  Point  Dynamic  Range  Waveform  for  m/  =  3,  sj 

1113  =  5,  S3  =  2. 
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Figure  4.4:  The  Point  Dynamic  Range  Waveform  for  mi  =  3,  s}  =  1,  m2  =  4,  s2  =  2, 

1713  =  5,  S3  =  0. 
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Figure  4.5:  The  Point  Dynamic  Range  Waveform  for  m,  =  3,  si  =  2, 1712  =  4,  S2  =  0, 

rri3  =  5,  S3  =  1. 
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Figure  4.6:  The  Point  Dynamic  Range  Waveform  for  m;  =  3,  sy  =  2,  m2  =  4,  s2  -  I, 

mj  =  5,  S3  =  0. 
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V.  APPLICATION  OF  THE  RSNS  TO  FOLDING  ANALOG-TO- 
DIGITAL  CONVERTERS 

A.         BACKGROUND  FOR  THE  FOLDING  ADC 

Analog  to  digital  converters  (ADCs)  are  critical  building  blocks  in  a  wide  range  of 
hardware,  from  radar  and  electronic  warfare  systems  to  multimedia  based  personal 
computer  and  work  stations.  The  best-known  architecture  for  high-speed  ADCs  is  the 
flash  converter  structure  [1.5].  In  this  structure  an  array  of  comparators  amplitude  analyze 
the  input  voltage  with  a  set  of  increasing  reference  voltages.  The  comparator  outputs 
represent  the  input  signal  in  a  thermometer  code,  which  can  be  easily  converted  into  a 
gray  or  binary  code.  This  architecture  requires  (2N  -  1)  comparators  to  achieve  TV-bit 
resolution,  which  makes  it  difficult  to  achieve  a  high  resolution  while  maintaining,  at  the 
same  time,  a  large  analog  bandwidth,  a  low  power  dissipation,  and  a  small  die  area. 

The  need  constantly  exists  for  converters  with  higher  resolution,  faster  conversion 
speeds,  and  lower  power  dissipation.  To  reduce  the  number  of  power  consuming 
components,  high  performance  ADCs  employ  a  parallel  configuration  of  analog  folding 
circuits  to  symmetrically  fold  the  input  signal  prior  to  quantization  by  high  speed 
comparators  (analog  preprocessing)  [6-9].  Recently  a  SNS  preprocessing  technique  has 
been  described  that  can  easily  be  incorporated  into  the  established  techniques  to  provide 
an  enhanced  resolution  capability  with  a  fewer  comparators  loaded  in  parallel.  The 
approach  is  based  on  preprocessing  the  analog  signal  with  symmetrical  number  systems 
(SNS  and  OSNS).  The  preprocessing  is  used  to  decompose  the  analog  amplitude  analyzer 
operation  into  /V  sub-operations  (channels)  which  are  of  smaller  computational 
complexity.  Each  sub-operation  symmetrically  folds  the  analog  signal  with  folding 
period  equal  to  the  channels  or  twice  the  channels.  Thus,  each  sub-operation  only 
requires  a  precision  in  accordance  with  that  modulus.  A  much  higher  resolution  is 
achieved  after  the  N  different  channels  are  used  and  the  results  of  these  low  precision  sub- 
operations  are  combined.  Increasing  the  resolution  of  these  architectures  can  be 
accomplished  by  incorporating  more  efficient  encoding.  One  of  the  important  factors 
which  helps  to  measure  the  efficiency  of  the  SNS  encoding  is  the  dynamic  range.    The 
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maximum  number  of  consecutive  vectors  without  a  redundancy.  The  other  factor  is  the 
total  number  of  comparators,  which  are  required  for  the  each  encoding.  The  7-bit  OSNS 
ADC  block  diagram  is  shown  in  Figure  5.1.  The  problem  with  the  SNS  and  OSNS  is  the 
encoding  errors.  The  transfer  function  of  7-bit  OSNS  and  the  encoding  errors  for  0% 
LSB  decimation  width  is  shown  in  Figure  5.2. 
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Figure  5.1:  The  7-bit  OSNS  Block  Diagram 
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Figure  5.2:  The  Transfer  Function  with  Isolation  Width  =  0%  LSB  for  7-bit  OSNS 
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B.         IMPLEMENTATION  ISSUES 

To  demonstrate  the  efficiency  of  the  RSNS  preprocessing  a  unipolar  5-bit  ADC 
with  mi  =  3,  m.2  -  4,  and  m3  =  5  is  considered.  The  dynamic  range  of  this  3-channel 
RSNS  is  43.  A  schematic  diagram  of  this  folding  ADC  is  shown  in  Figure  5.3. 

The  input  signal  is  folded  in  parallel  with  each  folding  period,  Prsns  equal  to  2Nm, 
where  N  is  the  number  of  channels  and  m  is  the  modulus  (channel  value).  The  output 
from  each  folding  circuit  is  quantitized  with  m  comparators.  For  this  simulation,  there 
will  be  3  +  4  +  5=  12  comparators  total  and  a  maximum  5  of  them  will  load  in  parallel. 
The  comparator  threshold  levels  are  adjusted  to  midlevel  quantize  the  folded  output 
waveform.  One  interconnecting  single  folding  and  comparator  stage  is  shown  in  Figure 
5.4  and  the  corresponding  values  of  the  resistors,  current  and  voltage  sources  are  shown 
in  Table  5.1.  The  encoder  to  a  binary  representation  for  a  high-resolution  digital  output 
converts  the  comparator  outputs  (the  RSNS  representation  of  the  input  signal). 

Each  RSNS  folding  circuit  must  symmetrically  fold  the  input  signal  at  least 

— ^-  number  of  times.  High  performance  folding  circuits  typically  consists  of  several 
2Nm 

identical  but  independent  stages  inter-connected  in  parallel  with  one  stage  for  each 

required  fold.  An  additional  stage  is  sometimes  used  to  register  the  input  voltages  in  DC. 

Vref  is  a  supplied  reference  voltage  to  each  stage  equal  to  the  input  voltage  point  at 

which  a  fold  (peak)  is  desired.    Also  shown  are  the  analog  input  node  and  the  folded 

output  node.   For  this  5-bit  example,  a  LSB  code  width  of  12mv  with  a  full-scale  voltage 

(26-l )  x  12mv  =  0. 756  v  is  chosen.  The  folding  period  for  the  modulus  m,  channel  is 

Tm.  =  2Nmi  x(LSB  Code  Width)  v  (53) 

where  N  is  the  number  of  modulus  and  m,  is  the  value  of  the  modulus  and  it  goes  from  1 
to  N.  The  first  required  reference  voltages  (yrefur  above  zero  and  Vref(or  below  zero)  are 

T 

Vrefo)  =  —  +  (s, x  (LSB  Code  Width))  v  (54) 
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Figure  5.3:  The  5-bit  RSNS  ADC  Block  Diagram 
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Resistor  Design 
Values 

RO 

noon 

Rl 

10,000  a 

R2 

900  Q. 

R3 

400  Q. 

R4 

300  Q 

R5 

300  Q. 

R6 

250  a 

R7 

250  Q  (  mod-3)  -  300  Q  (mod-4)  -  350  Q.  (mod-5) 

Voltage  Sources 
Design  Values 

VI 

15v 

V2 

15  v 

V3 

15v 

Current  Sources 
Design  Values 

11 

0.001  A 

12 

0.001  A 

Table  5.1:  The  Design  Values  of  the  5-bit  RSNS  (m;  =  3,  m2  =  4,  m3  =5) 


Vrefiv  = — —+(S:X(LSB  Code  Width))  v 


ref(O) 


(55) 
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where  Tm    is  the  folding  period  and  s,  is  the  shift  values  of  the  modulus  ra,.    The  other 

required  reference  voltages  are 

T 

Vref(N)  =  —  +  (s,x(LSB  Code  Width))  +  NTmi      v  (56) 

where  TV  is  the  N"7  reference  voltages  and  lies  between  {0,  1,2,  . . . }  for  the  modulus  ra,. 

The  folding  period  by  (53)  for  the  ra;  =  3  is  T3  =  2  x  3  x  3  x  72mv  =  276  rav.  That 

is,  the  folding  period  is  scaled  to  the  modulus.  For  ra/  =  3,  the  shift  value  is  taken  sj  =  0. 

T 
The  first  required  reference  voltage,  above  zero  is  then  Vre/(1)  =  —  =  108  rav  by  (54).  The 

reference  voltage,  below  zero  is  Vref{0)=-  —  =  -108mv  by  (55).    The  other  required 

reference  voltages  for  mi  =  3  are  calculated  by  (56)  and  shown  in  Table  5.2. 

The  folding  period  by  (53)  for  the  m2  =  4  is  T4  =  2  x  3  x  4  x  72mv  =  288  mv.  That 

is,  the  folding  period  is  scaled  to  the  modulus.  For  m2  =  4,  the  shift  value  is  taken  s2  =  1. 

The         first         required         reference         voltage,         above         zero         is         then 

T 
Vref(l)  =—+(slx(LSB  Code  Width))  =  \56mv  by  (54).     The  reference  voltage,  below 

T 

zero  is  Vref(0)  = — -  +  (s2x(LSB  Code  Width))  =  -\32mv  by  (55).    The  other  required 

reference  voltages  for  m2-4  are  calculated  by  (56)  and  shown  in  Table  5.2. 

The  folding  period  for  the  m3  =  5  is  T5  =  2  x  3  x  5  x  72mv  =  360  mv  by  (53).  That 
is,  the  folding  period  is  scaled  to  the  modulus.  For  m?  =  5,  the  shift  value  is  taken  S3  -  2. 
The  first  required  reference  voltage,  above  zero  is 

T 

thenVref(l)=  —  +  (s3x(LSB  Code  Width))  =  204mv    by   (54).      The   reference   voltage, 

T 
below  zero  is   Vref(0)= — -  +  (s3x(LSB  Code  Width))  =  -\56mv.     The  other  required 

reference  voltages  for  m3=  5  are  calculated  by  (56)  and  shown  in  Table  5.2. 
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Fold 
Number 

Modulus  3 
(T3=0.216) 

Modulus  4 
(T4=0.288) 

Modulus  5 
(T 5=0. 360) 

0 

-0.108 

-0.132 

-.0156 

1 

0.108 

0.156 

0.204 

2 

0.324 

0.444 

0.564 

3 

0.540 

0.732 

0.924 

4 

0.756 

1.020 

1.284 

5 

0.972 

1.308 

1.644 

6 

1.188 

1.596 

- 

7 

1.404 

- 

- 

8 

1.620 

Table  5.2:  The  Reference  Voltages  of  the  5-bit  RSNS  Design 

A  SPICE  simulation  is  used  to  generate  the  folding  waveforms  and  comparator 
outputs  states.  A  logic  map  (or  an  algorithm  (spice.m)  to  give  a  decimal  output  as  shown 
in  the  Appendix  processes  the  sampled  folding  waveforms  and  comparator  states.  The 
steady-state  folding  waveforms  for  each  channel  are  shown  in  Figures  5.5,  5.6,  and  5.7. 
These  folding  waveforms  are  generated  using  a  step  size  of  Av  =  0.00 1  v  in  the  SPICE 
simulation. 
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The  comparator  threshold  voltages  V,  for  each  modulus  derived  from  the  folding 
waveforms  to  mid-level  quantize  the  input  signal  into  RSNS  format.  The  threshold  levels 
occur  at  the  code  transition  points  and  are  tailored  to  the  particular  folding  waveform 
being  instrumented.  The  comparator  threshold  voltages  for  the  5-bit  RSNS  are  shown  in 
Table  5.3.  Note  the  threshold  values  are  not  uniformly  spaced.  The  comparator  outputs 
are  shown  in  Figures  5.5,  5.6  and  5.7  with  their  corresponding  folded  waveform  for  5-bit 
RSNS. 


Comparator 
Number 

Modulus-3 

Modulus-4 

Modulus-5 

1 

5.4097 

5.2994 

5.2234 

2 

5.6748 

5.5169 

5.4045 

3 

5.9452 

5.7408 

5.5946 

4 

- 

5.9668 

5.7878 

5 

- 

- 

5.9818 

Table  5.3:  The  Comparator  Threshold  Voltages  for  the  5-bit  RSNS  Design 

C.         TRANSFER  FUNCTION 

To  obtain  a  more  convenient  digital  output,  the  RSNS  encoded  signals  from  each 
channel  are  recombined.  The  comparator  outputs  from  each  channel  are  decoded  using  a 
logic  block  or  programmable  logic  array  (PLA).  The  RSNS-to-decimal  mapping  function 


50 


-0.2 


MOD3 


0.6 


20 


O       0 


-20  ■— 
-0.2 


-0.1 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


20 


-20 
-0.2 


-0.1 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


20 


O       0- 


-20 


-0.2  -0.1  0  0.1  0.2  0.3  0.4  0.5  0.6 

ADC  INPUT  VOLTAGE 

Figure  5.5:  Folding  Output  Waveform/Comparator  Outputs  for  the  RSNS  ray  =  3 
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Figure  5.7:  Folding  Output  Waveform/Comparator  Outputs  for  the  RSNS  m3  =  5. 
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for  the  5-bit  RSNS  design  is  shown  in  Table  5.4  where  each  entry  indicates  the  number  of 
comparators  in  the  ON  state. 


The  Point  Index 

RSNS  Moduli 

3 

4 

5 

0 

0 

0 

0 

1 

1 

0 

0 

2 

1 

1 

0 

3 

1 

1 

1 

4 

2 

1 

1 

5 

2 

2 

1 

6 

2 

2 

2 

7 

3 

2 

2 

8 

3 

3 

2 

9 

3 

3 

3 

10 

2 

3 

3 

11 

2 

4 

3 

12 

2 

4 

4 

13 

1 

4 

4 

14 

1 

3 

4 

15 

1 

3 

5 

16 

0 

3 

5 

17 

0 

2 

5 

18 

0 

2 

4 

19 

1 

2 

4 

20 

1 

1 

4 

Table  5.4:  The  5-bit  RSNS-to-Decimal  Mapping 
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The  RSNS  ADC  architecture  is  simulated  using  a  sampling  period  of  Av  =  0.001  v 
with  the  corresponding  transfer  function  shown  in  Figure  5.8.  The  LSB  code  width  is 
assumed  12  mv  for  this  simulation.  Therefore,  the  dynamic  range  lies  between  the  input 
voltages  of  0.012  x61  =  0.732  v  and  0.012  x  103  =  1.236  v.  The  entire  graph.  Figure  5.2 
lies  in  the  dynamic  range.  This  corresponds  to  the  portion  of  Figure  5.8  where  the 
dynamic  range  lies  from  0.732  v  to  1.236  v.  When  the  sampling  period  is  decreased  to  Av 
=  0.0004  v,  the  transfer  function  does  not  change.  Note  that  the  encoding  errors  typically 
present  in  the  SNS  and  OSNS  ADCs  do  not  exist.  This  clearly  demonstrates  the 
robustness  of  the  RSNS  encoding  (Gray  code  properties).  For  this  system,  channel  and 
shift  values,  the  dynamic  range  lies  between  the  point  indexes  61  and  103  as  shown  in 
Table  2.4. 

The  comparator  states  of  this  architecture  are  simulated  to  the  output  decimal 
value  by  using  the  algorithm,  trans.m.  The  beginning  point  index  (61)  of  the  dynamic 
range  is  assigned  the  decimal  value  0.  The  dynamic  range  point  indexes  are  assigned  the 
decimal  values  from  0  to  42.  The  point  indexes  within  the  dynamic  range  are  assigned 
the  decimal  values  starting  from  0  to  42.  Then,  starting  from  the  point  index  60  to  0,  the 
decimal  values  are  assigned  backwards.  If  the  point  index  of  the  discrete  states  is  a 
redundancy  within  the  dynamic  range,  it  is  assigned  the  corresponding  decimal  value.  If 
not,  it  receives  a  new  decimal  value  starting  from  -1.  Finally,  the  same  procedure  is 
applied  to  the  point  index  starting  from  104,  but  the  values  of  the  decimal  codes  start 
from  43.  The  dynamic  range  (useful  code  combinations)  and  the  redundancies 
(ambiguities)  can  be  seen  from  Figure  5.8. 
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Figure  5.8:  The  5-bit  RSNS  ADC  Transfer  Function  (Av  =  0.001  and  Av  =  0.0004  v) 
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VI.      CONCLUDING  REMARKS 

The  comparison  of  the  RSNS  with  the  other  SNS  formulations  is  shown  in  Table 
6. 1 .  The  folding  period  for  the  RSNS  is  a  factor  of  2N  greater  then  the  SNS  and  TV  greater 


Scheme 

Folding  Periods 

Number  of 
Comparators 
Per  Channel 

Dynamic  Range 

SNS 

m, 

_2_ 

N 

Y\mt 

SNS 

(all  moduli 

odd) 

m, 

_2_ 

min  j-n^,.  +-]!#!!,.[ 

SNS 

(one  moduli 

even) 

m, 

_2_ 

min  |~n^,  +  YlmA 

{   Z     1  =  2                l=j+\        J 

OSNS 

2ra, 

m, ;—  1 

N 

rim 

i=i 

RSNS 
(2  moduli) 

2Nmt 

mi 

4m,]  +  5ni2  -  5 
(moduli  1  or  2  apart) 

RSNS 
(2  moduli) 

2Nrrii 

m. 

4m\  +  5rri2  -  2 
(moduli  3  or  more  apart) 

RSNS 
(3  moduli) 

2Nm, 

mi 

3    2     15 

-m,  +—m,+7 

(moduli  1  apart) 

RSNS 

(3  moduli) 

2Nm, 

THi 

3    .     21 

— m,"  H rr\  4-14 

(moduli  2  apart) 

Table  6.1:  The  Summary  of  the  SNS  Systems 
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then  the  OSNS.  The  number  of  comparators  per  channel  for  the  RSNS  is  the  maximum  of 
all  the  SNS  systems.  The  dynamic  range  for  the  RSNS  is  less  than  the  OSNS,  but  greater 
than  the  SNS. 

The  essential  contribution  of  this  thesis  is  defining  a  new  symmetrical  number 
system,  robust  symmetrical  number  system  (RSNS)  and  investigating  its  properties.  High 
Speed  Analog-to-Digital  Converters  (ADCs)  for  2-channel  and  3-channel  are  summarized 
in  Tables  6.2  and  6.3,  respectively. 


Number 
System 

Channel 
(Modulus) 

Dynamic 
Range 

Folding 
Period 

Number  of 
Comparators 

SNS 

4 

1 

4 

2 

5 

5 

.2 

OSNS 

4 

20 

8 

3 

5 

10 

4 

RSNS 

4 

21 

16 

4 

5 

20 

5 

Table  6.2:  2-Channel  (modulus)  Summary  of  the  Number  Systems  for  mi  =  4  and 

m2  =  5 

The  dynamic  range  of  the  2-channel  RSNS  is  greater  than  the  SNS  and  OSNS,  but 
it  has  an  increased  number  of  comparators.  For  the  3-channel  RSNS,  although  the 
dynamic  range  is  greater  than  the  SNS,  it  is  less  than  the  OSNS.  The  RSNS  has  Gray  code 
properties  in  its  behavior.  Gray  code  properties  mean  only  one  comparator  threshold  is 
crossed  at  a  time  across  parallel  channels  (moduli).  This  property  of  the  RSNS  eliminates 
the  possibility  of  large  errors  at  the  code  transition  points  and  the  requirements  for 
interpolation  signal  processing. 

Using  SPICE,  explained  in  Chapter  V,  simulates  the  3-channel  RSNS  for  specific 
channel  and  shift  values.  From  the  simulation  results,  it  can  be  seen  that  the  RSNS  does 
not  have  any  encoding  errors.  This  property  of  the  RSNS  makes  itself  an  attractive  and 
useful  number  system  in  High  Speed  ADCs. 
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Number 
System 

Channel 

(Modulus) 

Dynamic 
Range 

Folding 
Period 

Number  of 
Comparators 

SNS 

3 

11 

3 

1 

4 

4 

2 

5 

5 

9 

OSNS 

3 

60 

6 

2 

4 

8 

3 

5 

10 

4 

RSNS 

3 

43 

18 

3 

4 

24 

4 

5 

30 

5 

Table  6.3:  3-Channel  (modulus)  Summary  of  the  Number  Systems  for  my  =  3,  m2  =  4 

and  rri3  =  5 

Future  efforts  will  attempt  to  resolve  the  following  important  questions. 

•  How  does  the  shift  effect  the  beginning  point  index  of  the  dynamic  range  for 
any  RSNS  configuration? 

•  Can  a  general  closed-form  expression  be  developed  for  the  3-channel  and  N- 
channel  RSNS? 
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APPENDIX 

A  .       DESCRIPTION  OF  THE  COMPUTER  ALGORITHM  DEVELOPED 

The  programming  in  this  study  is  accomplished  using  the  MATLAB  software 
package.  The  main  program,  which  is  used  to  find  the  properties  of  the  RSNS,  is  the 
rsns.m.  The  program  spice.m  and  trans.m  are  used  to  plot  the  SPICE  simulation. 

The  program  is  rsns.m  finds  and  displays  the  dynamic  range,  fundamental  period 
and  beginning-ending  point  index  of  the  discrete  states  without  redundancy  of  any 
number  of  RSNS  channels.  The  channel  and  shift  values  must  be  entered  to  run  the 
program.  This  program  also  finds  and  displays  the  point  dynamic  range  of  the  RSNS. 

The  program  spice.m  reads  the  folding  waveform  and  comparator  outputs  data  and 
displays  the  folding  waveform  and  comparator  outputs  for  each  channel  (modulus)  of  the 
RSNS  by  using  SPICE  simulation. 

The  program  trans.m  also  reads  the  folding  waveform  and  comparator  output  data 
which  are  calculated  by  using  SPICE  simulation.  The  algorithm,  trans.m  finds  the 
decimal  output  codes  for  each  comparator  outputs  of  the  RSNS  and  displays  the  transfer 
function  of  it.  The  dynamic  range  and  its  location  can  be  seen  in  this  algorithm  display 
(figure). 
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B.      PROGRAM  LISTING 

1.  RSNS  Properties  (rsns.m) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  rsns.m  % 

%  Modified  July  05th  1 997  % 

%  This  program  finds  the  dynamic  range,  point  dynamic  range,  fundamental  period     % 

%  and  beginning-ending  values  of  ambiguity  vectors  for  N-channel  RSNS  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

clear  all 

%  Enter  the  channel  number  and  the  values  of  them 

disp('This  program  finds  the  maximal  strings  of  non-redundant  vectors'); 

disp('for  the  N  channel  ROBUST  SNS'); 

%  Define  a  variable  for  the  channel  number 

chanum=input('Enter  the  Number  of  Channels  for  ROBUST  SNS  »  '); 

%  Find  the  fundamental  period  with  respect  to  the  entered  channel  and  shift  values 

period=  1 ; 

for  i=l:chanum 

m(i)=input(['Enter  '  int2str(i)  '.Channel  Value  »  ']); 

period=LCM(period,(2*chanum*m(i))); 
end 

nsearch=period+30; 
prompt='y'; 

%  Enter  the  corresponding  shift  values 
while  (prompt=='y')l(prompt=='Y') 
for  i=l:chanum 

s(i)=input(['Enter  '  int2str(i) '.  Channel  Shift  Value  »  ']); 
end 
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%  Initiate  the  variables  to  zero 
i=0; 
ii=0; 

j=0; 
jj=0; 

k=0; 

%  Define  a  variable  for  the  fundamental  period 

funper=0; 

%  Define  a  variable  for  the  dynamic  range 

dynrange=0; 

%  Define  the  waveform  for  the  entered  channel  and  shift  values  and  put  the  data  in 
%  matrix  "g" 
for  r=  1  :chanum 

%  Define  the  matrix  "mm"  to  put  the  entered  channel  and  shift  values 
mm(r,[12])=[m(r)s(r)]; 
%  Use  the  definition  of  RSNS 
for  i=l+s(r):chanum*m(r)+s(r) 
g(r,i)=floor((i-s(r))/chanum); 
end 
for  i=chanum*m(r)+l+s(r):  2*chanum*m(r)+s(r) 

g(r,i)=floor((2*chanum*m(r)+chanum-i+s(r)- 1  )/chanum); 
end 
%  Extent  enough  the  data  of  the  folding  waveforms  of  the  RSNS  to  find  the  dynamic 
%  range  and  fundamental  period 

g(r,2*chanum*m(r)+s(r)+ 1 :4*chanum*m(r)+s(r))=g(r,  1  +s(r):2*chanum*m(r)+s(r)); 
g(r,4*chanum*m(r)+s(r)+l:8*chanum*m(r)+s(r))=g(r,l+s(r):4*chanum*m(r)+s(r)); 
g(r,8*chanum*m(r)+s(r)+l:16*chanum*m(r)+s(r))=g(r,l+s(r):8*chanum*m(r)+s(r)); 
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g(r,  1 6*chanum*m(r)+s(r)+ 1 :32*chanum*m(r)+s(r))=g(r,  1  +s(r):  1 6*chanum*m(r)+s(r)); 
g(r32*chanum*m(r)+s(r)+l:64*chanum*m(r)4-s(r))=g(r,l+s(r):32*chanum*m(r)+s(r)); 
g(r,64*chanum*m(r)+s(r)+l:128*chanurn*rn(r)+s(r))=g(r,l+s(r):64*chanum*rn(r)+s(r)); 
end 

%  Define  a  new  matrix  "ga"  which  is  the  transpose  of  matrix  "g" 
ga=g'; 

%  Define  a  new  matrix  "gb"  which  is  extended  one  more  column  to  show  the  number  of 

%  rows 

gb(:,[2:(chanum+l)])=ga(:,[l:chanum]); 

[sgbr,sgbc]=size(gb); 

gb(:,l)=(l:l:sgbry; 

%  Define  a  new  matrix  "gc"  that  gives  the  row  number  through  the  number  of  search 

gc=gb; 

gc(nsearch+l:sgbr,  :)=[]; 

[sgcr,sgcc]=size(gc); 

%  Find  the  fundamental  period  to  check  this  result  with  the  formula  result  of  it 
xper=gc(  1  ,[2:(chanum+l )]); 
for  jj=chanum:sgcr; 

if  gc(jj,[2:(chanum+l)])==xper 
funper=gc((jj-(chanum- 1 )),  1 ); 
break 
end 
end 

%  Find  the  first  redundancies  in  matrix  "gc" 
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k=l; 

for  ii=chanum-l  msearch; 
xrec=gc(ii,[2:(chanum+l )]); 
for  ij=l  :floor((nsearch-ii)/(2*chanum)) 
jj=ii+2*chanum*ij; 

if  gc(jj,[2:(chanum+l)])==xrec 
redun=gc(jj,l); 

%  Define  a  new  matrix  "h"  to  put  the  first  redundancies 
h(k,l)=ii; 
h(k,2)=redun; 
k=k+l; 
break 
end 
end 
end 

%  Define  a  new  matrix  "  hsort"  to  sort  the  redundancy  matrix  "h" 

hsort=h; 

[yoy,ioi]=sort(hsort); 

%  Define  a  new  matrix  "  hsorted"  to  put  the  actual  row  number  of  "hsort" 

hsorted=[yoy(ioi(:,2),l)yoy(:,2)]; 

hsorted; 

%  Define  a  new  matrix  "hreduced"  to  eliminate  the  rows  of  the  matrix  "hsorted"  and  that 

%  will  not  allow  the  first  column  to  be  monotone  increasing 

[ssr,ssc]=size(hsorted); 

hreduced=hsorted; 

a=h(ssr,l); 
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[rx  cx]=size(hreduced); 
for  k=l:ssr 
for  i=l:ssr 
if  i<rx 
if  hreduced(i,l)==a 
hreduced(i+l  :ssr,  :)=[]; 
break 

elseif  hreduced(i+ 1 , 1  )<hreduced(i,  1 ) 
hreduced(i+l,  :)=[]; 
break 
end 
end 
[rx  cx]=size(hreduced); 
end 
end 
hreduced; 

%  Define  a  new  matrix  "H"  that  shows  the  point  index  of  the  distinct  vectors  without 
%  redundancy  and  their  lengths 
[hsr,hsc]=size(hreduced); 
H(l,l)=(chanum-1); 
H(2:hsr+l,l)=hreduced(l:hsr,l)+l; 
H(  1  :hsr,2)=hreduced(  1  :hsr,2)- 1 ; 
H(hsr+ 1 ,2)=nsearch; 

H(l:hsr+l,3)=H(l:hsr+l,2)-H(l:hsr+l,l)+l; 
[Hsr,Hsc]=size(H); 

%  Find  the  point  dynamic  range  of  the  RSNS  and  put  these  values  to  the  matrix  "HH" 
t=0; 
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for  n=l:Hsr 

forp=H(n,l):H(n+l,l) 
HH(p,l)=p; 
HH(p,2)=H(n,2); 
HH(p,3)=H(n,3)-t; 
t=t+l; 
ift=(H(n+l,l)-H(n,l)+l) 

t=0; 
end 
end 
n=n+l; 
if  n==Hsr 

break 
end 
end 

[HHsr,HHsc]=size(HH); 
HH(HHsr,:)=[]; 

%  Plot  the  point  dynamic  range  for  the  entered  channel  and  shift  values  of  RSNS 

figure(l) 

plot(HH(:,l),HH(:,3)),grid 

xlabel(THE  SHIFTED  CHANNEL  VALUES'),ylabel(THE  DYNAMIC  RANGE  OF 

THE  CHANNEL  VALUES') 

title([THE  DYNAMIC  RANGE  DISTRIBUTION  OF,num2str(chanum),'-CHANNEL 

RSNS']) 

orient  tall 

%print 

%  Find  the  dynamic  range  of  the  entered  channel  and  shift  values  of  RSNS 
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HHH=max(H); 
dynrange=HHH(:,3); 

%  Display  the  dynamic  range,  fundamental  period,  entered  channel  and  shift  values  and 

%  matrix  "H"  which  shows  the  %  beginning-ending  point  index  of  the  distinct  vectors 

%  without  redundancy  and  their  lengths 

dispC  ') 

disp([THE  FUNDEMANTAL  PERIOD  IS  ',num2str(funper),' ']) 

dispC  ') 

disp([THE  DYNAMIC  RANGE  IS  ',num2str(dynrange),' ']) 

fprintf('\n     BEGIN  POINT  INDEX  END  POINT  INDEX  DYNAMIC  RANGEVn) 

fprintf('%  1 1  .Of  %  1 6.0f  %  1 2.0f  \n',H') 

fprintf('\n     CHANNEL  VALUES      SHIFT  VALUESW) 

fprintf('%l  l.Of  %17.0f  \n',mm') 

prompt=input('Would  you  like  to  try  another  shift  (y/n)  ?  »','s'); 

end 

2.  Digital  Signal  Processing  for  the  RSNS  ADC  (spice.m) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  spice.m  % 

%  Modified  July  10th  1997  % 

%  This  program  finds  the  folded  waveform  and  comparator  outputs  of  % 

%  the  3-channel  RSNSwhich  is  simulated  in  SPICE  program  % 

%  mi  =  3         sj  =  0  % 

%  m2  =  4         s2  =  1  % 

%  m3  =  5         s3  =  2  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
clear  all 

%  Load  the  folding  circuit  and  comparator  outputs  of  the  modulus-3  RSNS 
load  newmod3.dat 
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a=newmod3; 

%  Plot  the  folding  waveform  and  comparator  outputs  for  the  modulus-3  RSNS 

figure(l) 

subplot(4, 1 , 1  ),plot(a(:,  1  ),a(:,2),,y*),grid,title(,MOD3,),ylabel(,FOLDED  OUTPUT') 

subplot(4,l,2),plot(a(:,l),a(:,5),'r'),grid,ylabel(,C3') 

subplot(4,l,3),plot(a(:,l),a(:,4),,g,),grid,ylabel('C2,) 

subplot(4, 1 ,4),plot(a(:,  1  ),a(:,3),'b'),grid,ylabel(,C  1  '),xlabel('ADC  INPUT  VOLTAGE') 

orient  tall 

%  Load  the  folding  circuit  and  comparator  outputs  of  the  modulus-4  RSNS 

load  newmod4.dat 

b=newmod4; 

%  Plot  the  folding  waveform  and  comparator  outputs  for  the  modulus-4  RSNS 
figure(2) 


subplot(5,l,l),plot(b( 
subplot(5,l,2),plot(b( 
subplot(5,l,3),plot(b( 
subplot(5,l,4),plot(b( 
subplot(5,l,5),plot(b( 
orient  tall 


,D,b( 
,l),b( 
,D,b( 
,l),b( 

,D,b( 


,2),'y'),grid,title(*MOD4'),ylabel('FOLDED  OUTPUT') 

,6),'m*),grid,ylabel('C4') 

,5),'r'),grid,ylabel('C3') 

,4),'g'),grid,ylabel('C2') 

,3),,b'),grid,ylabel('Cl'),xlabel('ADC  INPUT  VOLTAGE') 


%  Load  the  folding  circuit  and  comparator  outputs  of  the  modulus-5  RSNS 

load  newmod5.dat 

c=newmod5; 

%  Plot  the  folding  waveform  and  comparator  outputs  for  the  modulus-5  RSNS 
figure(3) 
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subplot(6, 1 , 1  ),plot(c(:,  1  ),c(:,2),'y' ),grid,title('MOD5,),ylabel(,FOLDED  OUTPUT') 

subplot(6, 1 ,2),plot(c(:,  1  ),c(:,7),,m,),grid,ylabel('C5') 

subplot(6, 1 ,3),plot(c(:,  1  ),c(:,6),,r,),grid,ylabel('C4') 

subplot(6,l,4),plot(c(:,l),c(:,5),'g,),grid,ylabel('C3') 

subplot(6,l,5),plot(c(:,l),c(:,4),'c'),grid,ylabel(,C2*) 

subplot(6, 1 ,6),plot(c(: ,  1  ),c(:,3),'b,),grid,ylabel('C  1  '),xlabel('ADC  INPUT  VOLTAGE') 

orient  tall 

%  Plot  the  folding  waveforms  for  modulus-3,  4,  and  5  RSNS  together 
figure(4) 

subplot(3,l,l),plot(a(:,l),a(:,2),y),grid,title('MOD3'),ylabel(,FOLDED  OUTPUT') 
subplot(3,l,2),plot(b(:,l),b(:,2),y),grid,title('MOD4'),ylabel('FOLDED  OUTPUT') 
subplot(3, 1 ,3),plot(c(:,  1  ),c(:,2),y),grid,title('MOD5'),ylabel('FOLDED  OUTPUT') 
orient  tall 

3.  Transfer  Function  for  the  RSNS  ADC  (trans.m) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  trans.m  % 

%  Modified  July  10th  1997  % 

%  This  program  finds  transfer  function  of  the  3-channel  RSNS,  which  is  simulated    % 
%  in    SPICE  program.  % 

%  m,  =  3         sj  =  0  % 

%  m2  =  4         s2=  1  % 

%  m3  =  5         s3  =  2  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
clear  all 
orient  tall 

%  Initiate  the  variables  to  zero 
i=0; 

j=0; 
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1=0; 

m=0; 

n=0; 

ka=l; 

kb=l: 

kc=l: 

kd=l: 

ke=l: 


%  Load  the  folding  circuit  and  comparator  outputs  of  the  modulus-3  RSNS  and  find 
%  howmany  comparators  are  on 
load  newmod3.dat 
a=newmod3; 
[sar,sac]=size(a); 
for  i=l:sar; 
al=a(i,3); 
a2=a(i,4); 
a3=a(i,5); 
if(al>0); 

A(ka,l)=l; 
else 

A(ka,l)=0; 
end 
if(  a2>0); 

A(ka,2)=l; 
else 

A(ka.2)=0; 
end 
if  (  a3  >  0); 
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A(ka,3)=l; 
else 

A(ka,3)=0; 
end 

ka=ka+l; 
end 
end 


%  Load  the  folding  circuit  and  comparator  outputs  of  the  modulus-4  RSNS  and  find 

%  howmany  comparators  are  on 

load  newmod4.dat 

b=newmod4; 

[sbr,sbcl=size(b); 

for  j=l:sbr; 

bl=b(j,3); 
b2=b(j,4); 
b3=b(j,5); 
b4=bG,6); 
if(bl>0); 

B(kb,l)=l; 
else 

B(kb,l)=0; 
end 
if  (  b2  >  0); 

B(kb,2)=l; 
else 

B(kb,2)=0; 
end 
if  ( b3>  0); 
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B(kb,3)=l; 
else 

B(kb,3)=0; 
end 
if(b4>0); 

B(kb,4)=l; 
else 

B(kb,4)=0; 
end 

kb=kb+l; 
end 
end 


%  Load  the  folding  circuit  and  comparator  outputs  of  the  modulus-5  RSNS  and  find 

%  howmany  comparators  are  on 

load  newmod5.dat 

c=newmod5; 

[scr,scc]=size(c); 

for  l=l:scr; 

cl=c(l,3); 
c2=c(l,4); 
c3=c(l,5); 
c4=c(l,6); 
c5=c(l,7); 

if(cl>0); 
C(kc,l)=l; 

else 
C(kc,l)=0; 

end 
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if  (c2>0); 

C(kc,2)=l; 
else 

C(kc,2)=0; 
end 
if(c3>0); 

C(kc,3)=l; 

else 

C(kc,3)=0; 
end 
if(c4>0); 

C(kc,4)=l; 
else 

C(kc,4)=0; 
end 
if(c5>0); 

C(kc,5)=l; 
else 

C(kc,5)=0; 
end 

kc=kc+l; 
end 
end 

%  Define  a  new  matrix  "D"  that  has  every  data  of  three  folding  circuit  and  shows  how 
%  many  comparators  are  on  or  off  in  each  of  them 
D(:,l)=a(:,l); 
for  m=l:sar; 

D(kd,2)=A(kd,  1  )+A(kd,2)+A(kd,3); 
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D(kd,3)=B(kd,  1  )+B(kd,2)+B(kd,3)+B(kd,4); 
D(kd,4)=C(kd,  1  )+C(kd,2)+C(kd,3)+C(kd,4)+C(kd,5); 
kd=kd+l; 
end 

%  Define  a  new  matrix  "DEC"  that  finds  the  transfer  function 
for  ke=l:sar; 
aa=D(ke,2); 

bb=D(ke,3); 
cc=D(ke,4); 
dd=D(ke,l); 
if((aa==0)&(bb==0)&(cc==  1 )) 

DEC(ke,l)=dd; 

DEC(ke,2)=-27; 
elseif((aa==O)&(bb==0)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-26; 
elseif  ((aa==l  )&(bb==0)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-25; 
elseif  ((aa==  1  )&(bb==  1  )&(cc=0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=33; 
elseif  ((aa==  1  )&(bb==  1  )&(cc==  1 )) 

DEC(ke,l)=dd; 

DEC(ke,2)=34; 
elseif  ((aa==2)&(bb==  1  )&(cc==  1 )) 

DEC(ke,l)=dd; 

DEC(ke,2)=35; 
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elseif  ((aa==2)&(bb==2)&(cc==  1 )) 

DEC(ke,l)=dd; 

DEC(ke,2)=6; 
elseif  ((aa==2)&(bb==2)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=7; 
elseif  ((aa==3)&(bb==2)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-24; 
elseif  ((aa==3)&(bb==3)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-23; 
elseif  ((aa==3)&(bb==3)&(cc==3)) 

DEC(ke,l)=dd; 

DEC(ke,2)=22; 
elseif  ((aa==2)&(bb==3)&(cc==3)) 

DEC(ke,l)=dd; 

DEC(ke,2)=23; 
elseif  ((aa==2)&(bb==4)&(cc==3)) 

DEC(ke,l)=dd; 

DEC(ke,2)=24; 
elseif  ((aa==2)&(bb==4)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-22; 
elseif  ((aa==l  )&(bb==4)&(cc=4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-21; 
elseif  ((aa==  1  )&(bb==3)&(cc==4)) 

DEC(ke,l)=dd; 
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DEC(ke,2)=-20; 
elseif  ((aa==l)&(bb==3)&(cc==5)) 


DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa==' 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa==' 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(kcJ 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 


)=dd; 
)=-19; 

0)&(bb==3)&(cc==5)) 
)=dd; 
)=-18; 

0)&(bb==2)&(cc==5)) 
)=dd; 
)=-17; 

0)&(bb==2)&(cc==4)) 
)=dd; 
)=-16; 

l)&(bb==2)&(cc==4)) 
)=dd; 
)=-15; 

l)&(bb==l)&(cc=4)) 
)=dd; 
)=15; 

l)&(bb==l)«fe(cc==3)) 
)=dd; 
)=10; 

2)&(bb==l)&(cc=3)) 
)=dd; 
)=41; 

2)&(bb==0)&(cc==3)) 

)=dd; 

)=-14; 


elseif  ((aa==2)&(bb==0)&(cc==2)) 
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DEC(ke,l)=dd; 

DEC(ke,2)=37; 
elseif  ((aa==3)&(bb==0)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=38; 
elseif  ((aa==3)&(bb==l)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=39; 
elseif  ((aa==3)&(bb==  1  )&(cc==  1 )) 

DEC(ke,l)=dd; 

DEC(ke,2)=-13; 
elseif  ((aa==2)&(bb==  1  )&(cc==  1 )) 

DEC(ke,l)=dd; 

DEC(ke,2)=35; 
elseif  ((aa==2)&(bb==2)&(cc==l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=6; 
elseif  ((aa==2)&(bb==2)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-12; 
elseif  ((aa=l  )&(bb=2)&(cc=0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=32; 
elseif  ((aa=l)&(bb==3)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-ll; 
elseif  ((aa==l)&(bb==3)&(cc==l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=28; 
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elseif  ((aa=0)&(bb=3)&(cc==l)) 


DEC(ke,l 
DEC(ke,2 

elseif  ((aa==< 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa==' 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 


DEC(ke,l 
DEC(ke,2 


=dd; 

=29; 

0)&(bb=4)&(cc=l)) 

=dd; 

=-10; 

0)&(bb==4)&(cc==2)) 

=dd; 

=-9; 

l)&(bb==4)&(cc==2)) 

=dd; 

=26; 

l)&(bb==3)&(cc==2)) 

=dd; 

=27; 

l)&(bb==3)&(cc==3)) 

=dd; 

=-8; 

2)&(bb==3)&(cc==3)) 

=dd; 

=23; 

2)&(bb==2)&(cc==3)) 

=dd; 

=42; 


elseif  ((aa==2)&(bb==2)&(cc==4)) 


=dd; 
=  19; 


elseif  ((aa==3)&(bb==2)&(cc==4)) 
DEC(ke,l)=dd; 
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DEC(ke,2)=20; 
elseif  ((aa=3  )&(bb==  1  )&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-7; 
elseif  ((aa==3)&(bb==  1  )&(cc=5)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-6; 
elseif  ((aa==2)&(bb==  1  )&(cc=5)) 

DEC(ke,l)=dd; 

DEC(ke,2)=17; 
elseif  ((aa==2)&(bb==0)&(cc==5)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-5; 
elseif  ((aa==2)&(bb==0)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-4; 
elseif  ((aa==l  )&(bb==0)&(cc=4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=14; 
elseif  ((aa==  1  )&(bb==  1  )&(cc=4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=15; 
elseif  ((aa=  1  )&(bb==  1  )&(cc=3)) 

DEC(ke,l)=dd; 

DEC(ke,2)=10; 
elseif  ((aa==0)&(bb==l  )&(cc==3)) 

DEC(ke,l)=dd; 

DEC(ke,2)=ll; 
elseif  ((aa==0)&(bb==2)&(cc==3)) 
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DEC(ke,l)=dd; 

DEC(ke,2)=-3; 
elseif  ((aa==0)&(bb==2)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-2; 
elseif  ((aa==  1  )&(bb==2)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=8; 
elseif  ((aa=l)&(bb==3)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=27; 
elseif  ((aa==l)&(bb==3)&(cc==l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-l; 
elseif  ((aa==2)&(bb==3)&(cc==l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=5; 
elseif  ((aa==2)&(bb=4)&(cc=l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=0; 
elseif  ((aa==2)&(bb==4)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=l; 
elseif  ((aa==3)&(bb==4)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=2; 
elseif  ((aa==3)&(bb==3)&(cc==0)) 

DEC(ke,l)=dd; 
DEC(ke,2)=3; 
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elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DECUe.l 
DEC(ke,2 

elseif  ((aa== 
DEQkeJ 


3)&(bb==3)&(cc=l)) 

=dd; 

=4; 

2)&(bb==3)&(cc==l)) 

=dd; 

=5; 

2)&(bb==2)&(cc==l)) 

=dd; 

=6; 

2)&(bb==2)&(cc==2)) 

=dd; 

=7; 


l)&(bb==2)&(cc==2)) 
=dd; 
=8; 

l)&(bb==l)&(cc==2)) 
=dd; 
=9; 

l)&(bb==l)&(cc==3)) 
=dd; 
=  10; 

0)&(bb==l)&(cc==3)) 
=dd; 
=  11; 

0)&(bb==0)&(cc==3)) 
=dd: 
=  12; 

=0)&(bb==0)&(cc==4)) 
)=dd; 
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DEC(ke,2)=13; 
elseif  ((aa==l)&(bb==0)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=14; 
elseif  ((aa==  1  )&(bb==  1  )&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=15; 
elseif  ((aa=  1  )&(bb==  1  )&(cc=5)) 

DEC(ke,l)=dd; 

DEC(ke,2)=16; 
elseif  ((aa=2)&(bb==l)&(cc==5)) 

DEC(ke,l)=dd; 

DEC(ke,2)=17; 
elseif  ((aa==2)&(bb==2)&(cc==5)) 

DEC(ke,l)=dd; 

DEC(ke,2)=18; 
elseif  ((aa==2)&(bb==2)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=19; 
elseif  ((aa==3)&(bb==2)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=20; 
elseif  ((aa==3)&(bb==3)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=21; 
elseif  ((aa==3)&(bb==3)&(cc==3)) 

DEC(ke,l)=dd; 

DEC(ke,2)=22; 
elseif  ((aa==2)&(bb==3)&(cc==3)) 
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DEC(ke,l)=dd; 

DEC(ke,2)=23; 
elseif  ((aa==2)&(bb==4)&(cc==3)) 

DEC(ke,l)=dd; 

DEC(ke,2)=24; 
elseif  ((aa==2)&(bb==4)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=25; 
elseif  ((aa==l)&(bb==4)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=26; 
elseif  ((aa==l)&(bb==3)&(cc=2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=27; 
elseif  ((aa==l)&(bb=3)&(cc==l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=28; 
elseif  ((aa=0)&(bb=3)&(cc=l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=29; 
elseif  ((aa=0)&(bb==2)&(cc=l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=30; 
elseif  ((aa==0)&(bb==2)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=31; 
elseif  ((aa=l  )&(bb==2)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=32; 
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elseif  ((aa==  1  )&(bb==  1  )&(cc=0)) 


DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 

elseif  ((aa== 
DEC(ke,l 
DEC(ke,2 


)=dd; 

)=33; 

l)&(bb==l)&(cc==l)) 
)=dd; 
)=34; 

=2)&(bb==l)&(cc==l)) 
)=dd; 
)=35; 

2)&(bb==0)&(cc==l)) 
)=dd; 
)=36; 

2)&(bb==0)&(cc==2)) 
)=dd; 
)=37; 

3)&(bb==0)&(cc==2)) 
)=dd; 
)=38; 

3)&(bb==l)&(cc==2)) 
)=dd; 
)=39; 

3)&(bb=l)&(cc==3)) 
)=dd; 
)=40; 

2)&(bb==l)&(cc==3)) 
)=dd; 
)=41; 


elseif  ((aa==2)&(bb==2)&(cc==3)) 
DEC(ke,l)=dd; 
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DEC(ke,2)=42; 
elseif  ((aa==2)&(bb==2)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=19; 
elseif  ((aa==  1  )&(bb==2)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-15; 
elseif  ((aa=  1  )&(bb==3)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-20; 
elseif  ((aa==  1  )&(bb==3)&(cc==5)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-19; 
elseif  ((aa==0)&(bb==3)&(cc==5)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-18; 
elseif  ((aa==0)&(bb==4)&(cc==5)) 

DEC(ke,l)=dd; 

DEC(ke,2)=43; 
elseif  ((aa==0)&(bb==4)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=44; 
elseif  ((aa==  1  )&(bb==4)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-21; 
elseif  ((aa==  1  )&Cbb==3)&(cc==4)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-20; 
elseif  ((aa==  1  )&(bb==3)&(cc==3)) 
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DEC(ke,l)=dd; 

DEC(ke,2)=-8; 
elseif  ((aa==2)&(bb==3)&(cc==3)) 

DEC(ke,l)=dd; 

DEC(ke,2)=23; 
elseif  ((aa==2)&(bb==2)&(cc==3)) 

DEC(ke,l)=dd; 

DEC(ke,2)=42; 
elseif  ((aa==2)&(bb==2)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=7; 
elseif  ((aa==3)&(bb==2)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-24; 
elseif  ((aa==3)&(bb==l)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=39; 
elseif  ((aa==3)&(bb==l)&(cc==l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-13; 
elseif  ((aa==2)&(bb==l)&(cc==l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=35; 
elseif  ((aa==2)&(bb==0)&(cc==l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=36; 
elseif  ((aa==2)&(bb==0)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=45; 
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elseif  ((aa==l  )&(bb==0)&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-25; 
elseif  ((aa==  1  )&(bb==  1  )&(cc==0)) 

DEC(ke,l)=dd; 

DEC(ke,2)=33; 
elseif  ((aa=l)&(bb==l)&(cc==l)) 

DEC(ke,l)=dd; 

DEC(ke,2)=34; 
elseif  ((aa==0)&(bb==  1  )&(cc=  1 )) 

DEC(ke,l)=dd; 

DEC(ke,2)=46; 
elseif  ((aa==0)&(bb==2)&(cc==  1 )) 

DEC(ke,l)=dd; 

DEC(ke,2)=30; 
elseif  ((aa==0)&(bb==2)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-2; 
elseif  ((aa==  1  )&(bb==2)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=8; 
elseif  ((aa==l)&(bb==3)&(cc==2)) 

DEC(ke,l)=dd; 

DEC(ke,2)=27; 
elseif  ((aa=l)&(bb=3)&(cc=3)) 

DEC(ke,l)=dd; 

DEC(ke,2)=-8; 
elseif  ((aa==2)&(bb==3)&(cc==3)) 

DEC(ke,l)=dd; 
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DEC(ke,2)=23; 
end 

ke=ke+l; 
end 
end 

%  Plot  the  transfer  function  of  the  3-channel  RSNS 

figure(l) 

plot(DEC(:,l),DEC(:,2)),grid 

xlabel('ADC  INPUT  VOLTAGE'),ylabel('ADC  DECIMAL  OUTPUT  CODE') 
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